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ABSTRACT 

This thesis investigates the estimation of Lanchester 
attrition- rat e coefficients in a deterministic aggregated 
combat model. Numerical values for the attrition-rate coef- 
ficients in the Lanchester-type model are taken to be the 
same as those for a corresponding continuous-time MarJcov- 
chain model, and maximum likelihood estimators are developed 
for these Markov-chain coefficients. A discussion of the 
historical background of Lanchester* s equations preceeds the 
theoretical development of components of the attrition rate 
equations. Using the functional forms and procedures devel- 
oped by Gordon M. Clark in his doctoral thesis, a simple 
computer simulation program is developed for a short dura- 
tion battle consisting of homogeneous forces. Recommenda- 
tions are made to modify this program to model heterogeneous 
forces in a similar battle. A discussion of the theoretical 
background supports these recommendations. 
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I* INTRODUCTION 



The purpose of this thesis is to investigate methods for 
estimating Lanchester attrition-rate coefficients, used in a 
deterministic aggregated combat model, from data generated 
by the detailed, high resolution Monte Carlo combat simula- 
tion model called STAR (Simulation of Tactical Alternative 
Responses) . An additional purpose is to construct a simple 
Lanchester-type combat model capable of using the estimated 
attrition rate coefficients to predict the outcome of STAR 
simulations. 

STAR is a model which has been developed primarily by 
O.S. Army students attending the Naval Postgraduate School. 
It has recently gained acceptance in the Army Community as a 
viable combat model, useful for studies such as the United 
States Army Mortar Study to be conducted in 1982. The cur- 
rent version of STAR is written in SIMSCRIPT and requires 
about 2 megabytes of core storage. The run time for one 
replication of STAR involving a company size force defending 
against a battalion size attacking force is approximately 
four CPU minutes on the IBM 3033 computer. 

The basic conceptual difficulties found in this thesis 
are summarized in Figure 1. The numerical values for the 
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attrition rate coefficients in the deterministic Lanchester- 



type model are assumed to be the same as those for a corre- 
sponding continuous-time Markov chain model. With this 
assumption, maximum likelihood estimators (MLE's) can be 
developed for the Markov chain model coefficients. Specific 
data needed for the HLE computations is extracted from the 
high resolution Monte Carlo simulation and used to compute 
•’estimated” coefficient values. These values for the 
continuous time Markov chain model are then substituted into 
the "analogous” deterministic Lanchester-type model. 

Constructing a Lanchester-type attrition model such as 
that proposed in this thesis is of significant value. With 
respect to the hierarchy of models concept [ Bef . 1], it 

could possibly fill the requirement as a ’hook* or link 
between STAB and low resolution models that exist at higher 
levels in the pyramid of models. However, possibly the most 
important contribution this model can make is as a filter 
model for STAB to reduce costs and time to run the many var- 
iations of a scenario in STAB. To explain further, consider 
the upcoming mortar study. To conduct such a study many 
runs/replications of a battle using the STAB model will have 
to be made. These numerous runs are required to compare 
various types of mortars, various types of force mixes. 
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Figure 1: Relationships of aodsls 

weapon sixes, and scenarios. Since STAR is a aonts Carlo 
siffluiation, var iance-reducing replications are required, 

1 1 




This requires significant computer time at a significant 
cost. By using a simpler analytical model which produces 
comparable output to a STAR run, many force mixes/weapon 
mixes/scenarios can be eliminated from further consideration 
at a minimal cost. STAR can then be used to concentrate 
only on the more promising options for force/weapons mixes. 

This thesis records the progress made in developing the 
small, simple, Lanchester- type combat model discussed above. 
It also lays the foundation for future work on this model. 
Chapter 2 of this thesis provides an introduction to Lan- 
chester-type combat models. Chapter 3 outlines the develop- 
ment of the Deterministic Aggregated Model of STAR on the 
Apple computer (DAHSTAC) , which represents the first step in 
building an analytical model comparable to STAR. Chapter 4 
explores various Lanchester attrition rate functional forms 
which may be used in the DAMSTAC or in future generation 
models. Chapter 5 details some desired enhancements of the 
model, while Chapter 6 discusses the work required to tran- 
sition from the current homogeneous force DAMSTAC model to a 
more realistic, albeit more complicated heterogeneous model. 
Chapter 7 contains the concluding remarks. 

Appendix A lists the variables used in the computer pro- 
gram for the DAMSTAC model. Appendix B lists some not.es of 
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interest concerning the program and model# which may not be 
apparent to the reader. Appendix C discusses the utility 
programs developed to support the model. Appendix D is a 
discussion of additional work done by the authors in the 
form of a computer model that explicitly plays the terrain 
effect. Appendix E consists of work done by two other stu- 
dents on a postprocessor which supports the computational 
requirements of the model and acts as an interface between 
STAR and the model. Finally, Appendix F contains the compu' 
ter listings of all programs written for this thesis. 
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II 



HISTORICAL BACKGROUND 



Lanchester-type combat models are analytical models in 
which differential equations are used to express the rate of 
change of force levels to represent the effects of the 
attrition process. The concept of representing attrition 
effects in land warfare with such mathematical equations was 
originated by F. W. Lanchester in 1914 in order to quantita- 
tively justify concentration of forces. 

Lanchester hypothesized that under conditions of "modern 
warfare" attrition between two homogeneous forces, each 
using "aimed" fire, could be modelled by: 



dx/dt = -ay (t) with x(0)= x<, 

dy/dt = -bx (t) with y(0)= y^ 

where x (t) and y (t) are the number of X and I firers alive 
at time t, and a and b are positive constants that are 
called Lanchester attrition- rate coefficients. The constant 
'a* represents the effectiveness of an x firer in killing y 
targets per unit of time and * b* represents the equivalent 
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effectiveness of the y firers. The equations written above 
are usually called "Lanchester ' s equations for modern war- 
fare”. Clearly, these functional forms for attrition rates 
that Lanchester developed are directly proportional to the 
attrition rate coefficients and the number of firers from 
the opposing force. The assumptions that must be made in 
order to use Lanchester *s aimed fire equations to model com- 
bat attrition are very restrictive. These assumptions can 
be stated as follows [Ref. 2]: 

1. Two homogeneous forces are engaged in combat; 

2. Each unit on either side is within weapon range of all 
units on the other side; 

3. The effects of successive rounds in the target area 
are independent; 

4. Each unit is sufficiently well aware of the location 
and condition of all enemy units so that it engages only 
live enemy units (one at a time) and hills them at a 
constant rate, which does not depend on the enemy force 
level. When an enemy target is hilled, search begins for 
a new target, with the rate of acquiring a new enemy tar- 
get being independent of the enemy force level; 

5. Fire is uniformly distributed over surviving enemy 
units. 
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Lanchester also hypothesized that when both sides used 
"area” fire, the attrition of forces could be modelled by: 

dx/dt = -axy 
dy/dt = -bxy 

The first three conditions under which Lanchester 's 
equations for area fire apply are identical to those for the 
equations for modern warfare, but assumptions 4 and 5 are 
replaced by [Ref. 3]; 

4. Each firing unit is aware only of the general area in 
which enemy forces are located and fires into this area 
without feedback about the consequences of its fire; 

5. Fire from surviving units is uniformly distributed 
over the area in which enemy forces are located; 

6. Each unit presents the same vulnerable area to enemy 
fire. This vulnerable area is much larger than the effec- 
tive area of a single round of enemy fire. Additionally, 
the number of hits required for a kill obeys a geometric 
probability law. 

The state equation that relates the force levels and 
does not explicitly contain time for the aimed-fire assump- 
tions is known as Lanchester 's Square Law, 
b (X2 - X2) = a (j;2 - yZ) . 
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The state equation for the area fire model yields 
Lanchester's Linear. Law, 

b (x^ - X) = a (Yo - y) . 

From these state equations many important results can be 
deduced, for example, battle outcome prediction conditions. 
Further analysis of Lanchester's classical equations and 
their implications will not be done here, suffice it to say 
that Lanchester's original work forms the basis for the ana- 
lytical combat model, with his simple model serving as a 
point of departure for the more complicated and elaborate 
analytical models discussed later in this thesis. Readers 
of this thesis and students who intend to continue with work 
initiated here should have a sound understanding of Lanches- 
ter's basic models before continuing further. 

Lanchester's equations for modern warfare and his equa- 
tions for area fire have, over the years, been thoroughly 
dissected and analyzed by students of differential equations 
and students of analytical combat models. One such student 
was Gordon Clark, a doctoral student at Ohio State Univer- 
sity in the late 1960 's. He was well aware of the power 
that simulation models such as Carmonette, a high resolution 
combat model, had for describing the complex process of com- 
bat. However, the major drawback to simulation, as he 
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perceived it, was that the method was slow and expensive. 
.Clark proposed that the use of an additional model to inter- 
pret simulation results would greatly benefit Army studies 
which required numerous runs with the simulation model. He 
went on to develop this analytical model which he called the 
Combat Analysis Model (COMAN) . This much simpler and faster 
running model would greatly reduce computer execution time 
and costs, and could serve as a filter model for Carmonette. 
A similar approach to Clark's methodology in developing 
COMAN will be used here to develop the analytical model for 
STAR. The next chapter details the first steps taken to 
build the analytical combat model for STAR. 
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III. DAaST AC MODEL 

The aggregated force-on- force analytical model is an 
integral part of combat modelling and serves several 
significant purposes. With the high cost of computer time 
and the limited computer core storage available to most mil- 
itary users, the analytical combat model offers a viable 
alternative to detailed high resolution simulations. 

Analytical combat models can be used in a stand alone 
configuration or in conjunction with a high resolution simu- 
lation. In the stand alone mode, the models are usually 
"transparent”, or, in other words, easy to understand. They 
are also small, guide but most importantly, they increase 
the user's understanding of the process of combat. In con- 
junction with high resolution combat simulations, the ana- 
lytical model offers, as Gordon Clark states in his doctoral 
thesis, the opportunity to "extrapolate the results of the 
combat simulation to predict the outcome from force mixes 
that were not explicitly simulated." In this way it can act 
as a filter to eliminate certain force mixes from further 
investigation. Also, Clark identifies another use for it as 
"an integral component of a large unit model to predict the 
outcome from individual unit actions." [Ref. U] 



19 



A major objective of this thesis was to build a simple 
aggregated combat model on the Apple II microcomputer. The 
input parameters for the model were estimated using the out- 
put data from a STAR run. This concept of estimating the 
attrition rate coefficients involved additional work in 
developing a data postprocessor compatible with STAR, and 
significant effort in the statistical analysis of the infor- 
mation provided by the postprocessor. 

There are several reasons why the Apple II computer was 
chosen for the development of the model rather than a main 
frame computer. First, the need exists for a 'portable' 
model that is easily accessible to the Army analyst. A desk 
top type system that can run a model quickly and with mini- 
mum setup time will be a valuable asset to the analyst. 

With the Apple II computer gaining wider acceptance in the 
Army for various uses, the accessibility of a model on such 
a system is greatly increased. 

Second, the nature of the model is such that it lends 
itself to a small system format. The aggregated form of the 
model reduces the required core storage from that required 
by a high resolution model. The concept of DAhSTAC is to 
serve as a filter model to weed out unacceptable force mixes 
or scenarios before turning to a high resolution model. The 
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idea here is to reduce the time and cost of running numerous 
iterations of a large, high resolution model on a main frame 
computer . 

Third, the cost of a computer system to support the 
model is certainly economical. At a cost of under $3000, 
the Apple II computer system with a model such as DAMSTAC 
provides a powerful analytical tool for a minimal cost, cer- 
tainly in relation to the cost of a main frame computer. 

Fourth, the Apple II has the capability to provide addi- 
tional analysis functions on the output of the model. With 
such programs as Apple Plot and available commercial statis- 
tical and numerical analysis packages, most of the needed 
analytical work can be accomplished on the Apple II compu- 
ter. This can result in additional savings in time by 
eliminating the requirement to transfer any data from one 
system to another. 

The laurels cast on the Apple II should not go without 
some discussion of its limitations. The Apple II is limited 
in capacity to 64 kilobytes of Random Access Memory storage 
space. So, essentially the model must be able to work 
within the amount of core storage. There are methods of 
working around this problem, if it arises, but this will 
usually create a new problem or limitation. For example. 
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some of the data or a part of the program can be kept on the 
storage diskette and put into con;puter memory only when ir 
is needed. This saves space in the computer’s memory but 
increases the run time of the model due to the time required 
to access the diskette to retrieve the required information. 

This raises another limitation of the Apple II — its run 
time. Compared to a main frame computer, the Apple II is 
very slow. This is, however, relatively speaking since com- 
puter time is normally measured in milliseconds and nanose- 
conds. For the analyst waiting for the output from DAMSTAC, 
the model on the Apple II would take about four to ten 
minutes to run. This is probably acceptable for most 
analysts. 

A. SCENARIO 

The model portrays a battle between a RED tank battalion 
attacking a BLUE tank company in the 10 x 10 KM Fulda Gap 
(Keyhole) terrain. See Figure 2 for the terrain map. 

The BLUE defensive positions and the RED positions along 
with their attack routes (coordinating points) are provided 
in Table 1. The RED units attack along three preplanned 
attack routes toward the BLUE defensive positions. Grid 
coordinates for unit locations and coordinating points were 
computed as the centroid of the elements of the unit. In 
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Figure 2: Terrain Map 
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Star's data base, each tank has its own individual location. 
For example, BLUE unit 1 has four different locations asso- 
ciated with the it, each one corresponding to a tank in the 
unit. The centroid of the unit was taken to be the average 
of these grid coordinates. Thus, the initial location for 
BLUE unit 1 in Table 1 is the average of these coordinates; 

54150 100690 54150 100620 

54190 100520 54160 100790 

The rate of RED advance is a constant but specified by the 
user. BLUE units can move to alternate positions at a rate 
designated by the user. However, for this scenario, BLUE 
units remain stationary for the entire battle. 

To obtain the raw output data that would correctly 
reflect the above scenario, it was necessary to ’’turn off*' 
many of the submodels of STAS. For example, the Field 
Artillery module was shut off by scheduling the first fire 
mission well after the battle ended. Weapon systems other 
than tanks within the fighting units were rendered useless 
by setting their ammunition counts (basic loads) to zero. 
Additionally, the communications and the Air Defense modules 
were not played. After all of these preliminary actions 
were taken, the result was a battle exclusively between 
tanks. This effort worked well to reduce a complicated. 
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heterogeneous force battle to a simple homogeneous force 
battle. One minor problem arose, however, when a weapon 
system other than a tank was detected and fired on when, by 
the scenario, it did not even exist. The effects of this 
target acquisition and engagement were considered insignifi- 
cant and were disregarded. 



Table 1: UNIT LOCATIONS/COORDINATING POINTS 







X-COORD 


Y-COORD 


SIZE 


REMARKS 


UNIT 


1 


54163 


100655 


4 


initial location 


UNIT 


2 


54788 


98120 


4 


initial location 


UNIT 


3 


57448 


97893 


4 


initial location 


BED 




X-COORD 


Y-COORD 


SIZE 


REMARKS 


UNIT 


1 


58727 

52000 


101491 

96000 


10 


initial location 
coordinating pt 


UNIT 


2 


581 10 
52000 


1016 96 
96000 


10 


initial location 
coordinating pt 


UNIT 


3 


58727 

52000 


101491 

96000 


10 


initial location 
coordinating pr 



B. ASSUMPTIONS 

The attrition rate functional form chosen for the 
DAMSTAC model, which will be discussed later in this thesis, 
requires the following assumptions. 

1. The conditional kill rates, a and 0 , are constant and 
independent of time if an individual firer has at least 
one acquired target. If there are no acquired targets. 
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ths kill rates are zero. 



2. P and q, the probabilities of nondetection, are 
assumed to be constant functions of time, that they are 
independent of the force size and the number of acquired 
targets, and that the probabilities p and q are equal for 
all x-firer/y-target and y-f irer/x-tar get combinations 
respectively. 

3. Synergistic effects are neglected, therefore, the 
actual force attrition rate is equal to the sum of the 
individual kill rates for the opposing forces. 

additionally, the assumptions listed below are designed to 
remove the distractions of a complicated scenario and to 
make the initial model easier to understand: 

1. All units have point locations; ie. a battalion is 
represented by a grid coordinate and all calculations for 
the battalion are made from that location. 

2. Units are homogeneous. That is, they are composed of 
identical firing elements/weapons systems. 

3. Fire is uniformly distributed over surviving enemy tar- 
gets within range. 

4. Sufficient supplies of ammunition for both sides exist 
for the duration of the battle. 

5. Units move at a constant rate of speed. 
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C. SYNOPSIS OF THE COMPUTER MODEL 



The DAMSTAC computer program is a modular type program 
written in the Applesoft BASIC language. The program 
accesses external data files as a source of input of initial 
data for the two forces. These data files are produced by 
the user with the help of the Force Maker utility program 
described in Appendix C. Additionally, the program has the 
capability to make output files that can be stored on a 
diskette and accessed by the user with the Results Header 
utility program, also described in Appendix c. Results 
Reader will print the output to the screen or to a printer. 

A copy of the program for the model is included as List- 
ing A in Appendix F. A sample run of the model including 
user inputs is included as Figure 3. Flow charts of the 
main program and subroutines are included as Figure 4, Fig- 
ure 6, and Figure 7. 

"I • Initializatio n of Arr ay s and V aria bl e s 

Lines 30 thru 170 dimension the arrays and matrices 
and initialize hardwired data. Dimensions are set to handle 
up to five units in each force and up to ten coordinating 
points, including the initial location, for each unit. The 
user can modify these capacities by changing the dimensions 
of the appropriate variables. 
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Lines 180 thru 350 read BLUE force data from a file 



to include the number of units, initial locations, any coor- 
dinating points, and rates of advance, and compute total 
force level based on user input of the unit force levels. 

The time increment (DT) is set to 30 seconds while the maxi- 
mum time of the battle (MT) is set to 1500 seconds. These 
values are, of course, adjustable by changing them in line 
80. For both forces, the breakpoint criteria is set at 0.3 
of the force level and movement criteria is set at 0.5 of 
the unit strength, even though the movement criteria was not 
used for any runs of the model. 

Lines 360 thru 520 perform the same operations for 
the P.ED forces as lines 180 thru 350 perform for the BLUE 
forces. 

Lines 530 and 600 ask if the printer should be 
turned on or off and calls the printout subroutine to print 
the initial data. 

2. Mis P rogram (lin es 61 0-78 0) 

The main program performs these operations: 
increments the time, 

checks if the time exceeds the max time, 

calls the movement and attrition subroutines, 

recomputes current total force levels after attrition, and 
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Figure 4; Flow Char"t of tiie Maiu Program 
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checks force levels against the designated breakpoints for 
each force, to determine if the battle has ended. 

If the breakpoints or the max time are exceeded, the battle 
terminates and the program prints the force levels. Addi- 
tionally the program will write the force levels for each 
time interval onto the diskette, (lines 1470 thru 1580) 

3, ao vemen t Subroutine ( lines 790- 1 160 ) 

Both RED and BLOE forces are capable of moving, 
aovement of the BLUE force is nullified by setting its move- 
ment rate to zero in the initialization part of the program. 
The movement subroutine uses basic geometry to compute new 
coordinates for the moving unit (see Figure 5) . The move- 
ment process is based on moving a unit from one coordinating 
point to the next in a straight line. The amount of move- 
ment depends on the rate of movement (BE (I) or RH (I) ) and 




Figure 5; Geometric Representation of Unit aovement 



31 




Figure 5; Plow Chart of the aoveoent Subroutine 
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the time interval length (DT) . If the distance computed is 
farther than the next coordinating point, only movement to 
the coordinating point is accomplished. Lines 800 to 970 
perform computations to move RED forces while Lines 980 to 
1150 perform the same computations for BLUE forces. 

'*• Attrit i on S ubrou tine ( line s 1 1?0-1480> 

A range check is made for each firer/target combina- 
tion to insure the units are within firing range. If so, 
attrition caused by that firer on the target is computed. 
Next total attrition for each unit is computed based on the 
partial attritions of the firers as a fraction of their 
force levels. In other words, when a firer is firing at 3 
enemy units, each target unit receives 1/3 of the firer's 
fire power. 

5- P rint out Subr o ut ine (lines 1490-1530 & li nes 
17 60-1 940) 

There are two different output formats used in the 
program. The first (lines 1490 to 1530) prints only a sum- 
mary of the total force levels for both forces for each time 
interval. This is the same data that is saved to a file by 
the final portion of the program. The second format (lines 
1760 to 1940) prints out the units of each force with their 
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Figure 7: Flow chart of the Attrition Subroutine 
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current locations and strengths. This format allows the 
user to see which unit, specifically, is being attritted and 
the grid coordinates of that unit at the time of attrition. 

6. Output File Making S ubroutine f lines 1540-1750 ) 

This subroutine allows the user to make a file of 
the output from the current run of the program. This 
routine actually produces two files, one consisting of the 
attrition rate coefficients and probabilities of non-acqui- 
sition, and the second containing the time and force levels 
of both forces at that time. These files are separated so 
that the second file can be used with the APPLEPLOT program 
to produce a graph of the two force levels as a function of 
time. This file has the postscript, ‘'.RESULTS”, added to 
its file name, while the first file with the attrition rate 
coefficients has the postscript ".COEFS". 

D. METHODOLOGY FOR COMPUTING PARAMETER ESTIMATES 

The process of building a computer simulation combat 
model requires the development of methodologies for 
modelling the various aspects of combat. Constructing a 
force-on-force attrition model requires methods for consoli- 
dating combat phenomena into simplified algorithms with 
coefficients representing one or more facets of combat. 

This section discusses one significant part of the 
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force-on-f orce attrition model: the methodology for estimat- 
ing the attrition rate coefficients in the Lanchester type 
attrition rate equations. 

In deriving the equations for the attrition rates (dx/dt 
and dy/dt ) , it is necessary to understand the processes that 
must occur to cause a firer to inflict attrition. These key 
processes can be summarized as: 
the environmental process; 
the target acquisition process; 
the target engagement process. 

The environmental process involves the terrain and weather 
factors that affect firer/target line of sight (LOS) . The 
target acquisition process represents the visible capability 
of the firer to detect a target on the battlefield, while 
the target engagement process describes rhe firer *s ability 
to hit and kill an acquired target. These three processes 
will all be represented by the attrition rate coefficient, 

4, in the attrition rate equation. This coefficient will be 
broken down into its three components which describe mathe- 
matically the phenomena of the three processes. The three 
components are: 

= the probability that a target is visible to the firer 
(environmental process) 
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= the probability that a visible target is acquired 
(target acquisition process) 

Q!= the casualty rate at which a firer attrits the 
acquired target (target engagement process) 

a should not be confused with A, the attrition rate 
coefficient, of which a is a factor. The attrition rate 
equation may be represented in words by: 

dx/dt = - [P(tgt is visible)] (P (tgt is acquired)] [casualty rate] 
The second probability above can be broken down further to 
P(tgt is acquired j tgt is visible). Note that this is a 
conditional probability since it is assumed that the target 
must be visible to be acquired. By identifying the 
mathematical expressions for each of these components of the 
attrition rate coefficient, the attrition rate equation can 
be obtained. 

I* The Bny ironmenta l P roc ess 

There are generally two ways to represent line of 
sight in a model: 

as a mathematical simulation of actual terrain; 
as a stochastic process. 

The mathematical simulation method normally involves a 
digital representation of terrain elevation over a segment 
of the terrain. The elevations of the segments between the 
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firer and the target are used to determine if line of sight 
exists. Virtually all high resolution models use this 
method, or variations of it, to determine line of sight. 

STAR uses a different form of mathematical simulation to 
represent terrain. This is called functional terrain, in 
which hills and forests are represented by elliptical con- 
tours of elevation. Hills are modelled by describing them 
in terms of bivariate normal equations with varying parame- 
ters for each hill. Forests or built-up areas are repre- 
sented by an elliptical pattern on the ground again with 
varying parameters for each forest. 

For a low resolution or highly aggregated model, rhe 
mathematical simulation method for line of sight is usually 
not practical or desireable in terms of core space and com- 
putation time required. Therefore, the second method, in 
which LOS is represented as a stochastic process, is used. 
The concept of this process is as follows. The effect of 
terrain on rhe process of line of sight is that a target 
falls into one of two categories or states: either it is 

visible or it is invisible. 

Assuming a steady-state situation for the terrain, 
probabilities can be assigned to each of these two states. 
These probabilities represent the percentage of time that 
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the target will be in that state (visible or invisible) for 
that terrain over a reasonable amount of time. Graphically 
the states of the target can be represented by the diagram 
in Figure 8. 




Since this is a stochastic process, there is an 
exponential distribution associated with each of the two 
states. This distribution represents the time between two 
occurrences of the same state. On the graph in Figure 8, 
the time between the end of t'| and the start of is expo- 
nentially distributed with mean, rj . Likewise, the time bet 
ween the end of and the start of is exponentially 
distributed with mean, n . Graphically, this can be repre- 
sented by Figure 9, where the circles represent the two 
states, visibility and invisibility, and the arrows 
represent the transition from one state to the other at the 
given rata. 
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Figure 9: Two-srate Stochastic Model of Intervisibility 

The attrition caused by a firer on a target can 
occur only when the target is visible (assuming no atxrition 
from wild, blind shots or indirect fire) . Therefore, the 
attrition rate is a function of the probability that the 
target is visible to the firer. This probability will be 
denoted as , and can be expressed as the proportion of 
time that the target is in the visible state in a steady 
state situation, p = .. • • 

2. The Target Acquisition Process 

Assuming now that line of sight exists between a 
firer and a target, the process of target acquisition must 
now be represented. There are two primary modes of target 
acquisition: 

Serial acquisition, and 
Parallel acquisition. 

In the serial acquisition process, targets are acquired 
between engagements. That is, in concept, the firer is 
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concentrating on the engaged target and does not or is not 
able to acquire other targets on the battlefield. Any addi- 
tional acquisitions made prior to engaging the current tar- 
get are lost and must be re-acquired after the current 
engagement has ended. Examples of firers who use serial 
acquisition are: 

any firer who is “buttoned up“ and has his field of view 

narrowed to a thin band; 

wire guided missiles/anti-tank systems. 

In the parallel acquisition mode, the firer main- 
tains his list of acquired targets throughout an engagement 
(except those lost through loss of line of sight) . Thus, 
when an engagement has ended, by the target being killed or 
loss of line of sight, the firer may immediately engage a 
target that was previously acquired and is still a valid 
target. All systems which do not use serial acquisition are 
assumed to use the parallel mode. 

For purposes of this thesis, all firers are assumed 
to use parallel acquisition. Therefore, serial acquisition 
will not be discussed any further. The reader should note 
that the serial acquisition process does have validity in 
modelling action on the battlefield, and should be consid- 
ered in future enhancements of this model. 
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To explain the derivation of the target acquisition 
portion of the attrition rate equation, the diagram in Fig- 
ure 9 must be modified to that in Figure 10. 



INyiSl8L£ 
STA7£ 

Figure 10: Three-state Stochastic Model of Intervisibiliry 

The two-state stochastic process has been changed no 
a three-state process by dividing the visible state into the 
visible-not acquired state (VHA) and the visible- acquired 
state (Vi) . What has occurred is that a target can move 
from the visible-unacguired state to the visible acquired 
state at a rate \ . Since both of these states comprise the 
visibility state in total, the rate of moving from either of 
them to the invisibility state remains at . Notice that 
there is no direct transition from the invisible state to 
the visible-acquired stats. This is only a matter of 
representation, since it is assumed that no more than one 
event, ie. moving to a new stats, will occur during a time 
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increment. Thus, the situation described can still be 
represented by Figure 10. 

Another point of interest in this diagram is the 
fact that there is no transition from visible-acquired to 
visible-not acquired. This implies the assumption that once 
a target is acquired, it remains acquired unless it is 
killed or line of sight is lost, in which case it transfers 
to the invisible state. 

To determine the effect of the target acquisition 
process on the the attrition rate, the probability that a 
target is acquired must be determined. From Figure 10, and 
using the same concept of a stochastic process as that used 
in section III-D-1, the equation can be written 
(tgt acquired) = ■ ^ ^ ^ . (1) 

This should more properly be stated as the probability that 
a target is acquired given that it is visible, since this is 
actually a conditional probability with the state of acqui- 
sition being dependent on the target being visible. 

Therefore, the probability that a target is visible 
and acquired in the steady state is: 



va 



P (tgt vis & acq ) 

P^ (tgt acq I tgt vis) P„(tgt 

a V 



(\+m) 



vis) 



( 2 ) 
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The implications of this equation will be discussed later in 
Chapter V. 

3* Formula t ing the A ttrition- rate Equation 

Using equation (2) in the attrition- rate equation, 
the attrition rate of a single firer on a single target can 
be summarized using expected values as: 

dx/dt = Of (3) 

That is, the attrition rate equals the rate of killing an 
acquired target times the probability that a target is visi- 
ble and acquired. However, since there are likely to be 
more than one firer and more than one target on the battle- 
field, this equation must be modified to reflect this. 

If a firer is going to cause attrition on his tar- 
gets, he must be able to acquire one or more of these tar- 
gets. For firer y1 with x targets, the probability that y1 
acquires one or more targets is just one minus the probabil- 
ity that he acquires no targets. Assuming the same 
probability of acquiring any of the targets, the probability 
of acquiring none of the targets is (1 - . Therefore, 

V a 

'■a,, = ' - n - • C) 
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The attrition rate now becomes 



dx/dt = -P a (5) 

for a single firer against x targets. Since there are more 
than one firer on the battlefield, each assumed to be acting 
independent in this process, the attrition rate in equation 
(5) must be multiplied by the number of firers(y) and the 
equation is now complete; 

dx/dt = -P Q! y . (6) 

In terms of the probability of acquiring a target, this 
equation can be written as; 

dx/dt = -(1 - (1 - p^^)^ ) Qty . (7) 

This is the same equation which Clark uses in his COMAN 
model, that uses results from a high fesolurion model to 
estimate the value for the conditional casualty rate and the 
probability that a target is not acquired. Henceforth, (1 - 
P^^) will be referred to as p for RED firer/BLUE target com- 
binations (or q for a BLUE firer/EED target combination) . 

4. The Target En g ageme nt P roc ess (Bonder vs Clark ) 

The target engagement process can be modelled in its 
most detailed form, as in the high resolution model, or in 
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its simplest form, as in the aggregated, analytical model. 
The high resolution model, such as STAR, includes explicit 
representation of every factor that is involved in the pro- 
cess of aiming at, firing at, hitting, and killing a target. 
These factors include: 

weapon system characteristics; 
weapon/ammunition type characteristics ; 

target characteristics such as size, shape, vulnerability, 
etc. 

In contrast, the low resolution, aggregated model 
represents the engagement process by one value — the attri- 
tion rate coefficient. This is the method to be employed in 
the model presented here. The task now is to develop a 
method to accurately estimate this coefficient for use in 
the attrition rate equation. Note that we are actually 
talking about two attrition rate coefficients and two equa- 
tions. The derivation and method for estimating one is the 
same as the other, so for purposes of the discussion, refer- 
ences to an attrition rate coefficient will be in the 
singular. 

There are presently two basic approaches to the 
problem of estimating attrition rate coefficients for paral- 
lel acquisition and stochastic line of sight. The first is 
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the use of the independent analytical model, involving 
weapon system characteristics and the use of algorithms to 
determine the expected time for a fireir to kill a single 
target. The attrition coefficients are then the reciprocal 
of these expected times. This method will be referred to as 
Bonder's method after Seth Bonder who has used it exten- 
sively in his combat models such as VECTOR II. 

The second method is the fitted- parameter analytical 
model, involving the use of a separate high resolution, 

Monte Carlo combat simulation to produce data to determine a 
maximum likelihood estimate of the attrition coefficients. 
This method will be referred to as Clark's method, since 
Gordon Clark popularized it in his doctoral thesis and in 
the development of the Combat Analysis Model (COMAN) . 

Most of this section involves a discussion of 
Clark's method and the practical application of it. How- 
ever, in order to provide a contrast to that method, a short 
description of Bonder's method will be provided. 

In general, the concept of Bonder's method is to 
produce an algorithm representing the expected time that it 
takes a firer to kill a single target. The elements of this 
algorithm include computations for line of sight, target 
acquisition, target selection, and the killing process. 
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Each alement accounts for a certain amount of time in the 
expected time to kill a target. This entire process is seen 
by Bonder as a Markov renewal process where the outcome of a 
round fired is dependent only on the outcome of the previous 
round. By using weapon systems characteristics data as raw 
input, times to acquire, select, and kill a target can be 
determined and used in the algorithm for expected time to 
kill a target. 

Bonder's method uses an independent, analytical 
model requiring a large data base of factors, such as range 
dependent weapon system characteristics, quantifications of 
line of sight and target acquisition parameters, to support 
the computation of the attrition rate coefficient. 

A detailed explanation of Bonder's method for esti- 
mating the attrition rate coefficients can be found in the 
basic documentation for the VECTOR II model [Ref. 5]. Basi- 
cally, Bonder says that the attrition rate coefficient is 

represented in the parallel acquisition mode by A^^j = at 

where Qj^jis the probability that a given group-i weapon is 
firing at a group-j target. Q_ can be expressed in terms 
of the probability that a group-j target is visible and has 
been detected by a grcup-i firer. This probabiltiy is 
denoted as S-. , and the expression for Q.. is : 
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Q = 5 /n (1 - s. . ) for j > 2 . 

ID ID ■'■J 

Assuming homogeneous forces, this equation reduces to 

Q.. = S. . , and the i-j subscripts can be eliminated 
ID ID 

since there is no distinction among types of force elements. 
His resulting expression for S is identical to that of equa- 
tion 4 , where S..and P_ are defined the same. 



ID 






Bonder's expression for the conditional casualty 
rate, a; » is written as the inverse of the expected time to 
kill an acquired target. This expression, when further bro- 
ken down contains elements dependent on the weapon system 
characteristics, target characteristics, and a number of 
other subsystem performance parameters. These parameters 
were determined empirically outside the operational environ- 
ment of the battlefield and molded together in this pre- 
determined structure zo produce the value of the attrition 
coefficient. This derivation of the attrition coefficient 
is the point where Bonder's method and Clark's method 
differ. 



Clark's methodology is based on the assumption that 
the process of attrition in combat is a nonhomogeneous, two 
independent type Poisson Process where the time between 
casualties has an exponential distribution with parameter \ , 
called the total conditional casualty rate. Likewise, the 
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time between RED casualties and the time between BLOE 
casualties each has an exponential distribution with parame 
ters Xj. and respectively , where \ = . Prom the 

theory of Poisson Processes, the probability of a RED casu- 
alty occurring is^r/;y^ , and the probability of a BLUE casu- 
alty occurring is^b/^^ . With this basis, Clark is able to 
develop a likelihood functon in terms of the conditional 
casualty rates a and /3 , of which the casualty rates are 
functions, ie«, Xj^= f (or) and Xj-= f (/J) • [Ref. 6] It should 
be noted that, in order to maintain consistency of terminol 
ogy, the equations attributable to Clark's work have been 
adjusted to reflect the definitions presented in this the- 
sis. 

Clark also uses the concept that attrition is a 
function of target acquisition, which depends on the prob- 
ability of a firer acquiring a target. This concept pro- 
duces two more parameters which must be estimated: p, the 
probability that a BLUE target is not acquired by a RED 
firer; q, the probability that a RED target is not acquired 
by a BLUE firer. 

For the functional form in Clark's COWAN model, the 
equations can be written as: 
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( 8 ) 



-Q![l-p"*] with m(o) = 
ff= -^(1- q" ] with n(o)= no 

where the parameters are defined by: 

fj^= the conditional attrition rate function for BLUE 
forces 

ff= the conditional attrition rate function for RED forces 
a- conditional kill rate for a RED weapon against acquired 
BLUE targets 

BLUE weapon conditional kill rate corresponding to a 

p = probability that a BLUE target is unacquired by an 
individual RED fixer 

q = probability that a RED target is unacquired by an 
individual BLUE firer 

m = number of BLUE fixers surviving 

n = number of RED firers surviving 

These functional forms depend directly on the abil- 
ity of firers to detect targets. To see this, consider just 
the conditional attrition rate function for the BLUE force, 
f^, . The definition of p implies that the value of {l-p"*) is 
the probability that a RED firer has at least one acquired 
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BLOE target. Therefore, the attrition rate for the BLUE 
force is the kill rate for acquired targets times the prob- 
ability of having an acquired target times the number of 
firers. 

It should be noted that Lanchester*s classical 
attrition rate functional forms are special cases of Clark's 
attrition rate functional forms, and, under appropriate con- 
ditions, Clark's equations reduce to Lanchester's equations. 
When p and q are zero, or very nearly zero, 

- oc (l-o"*) n = - a n (9) 



and. 



f^ = - 0 ( 1-0'') m = - 0 m (10) 

which are Lanchester's classical aimed fire equations. When 
individual targets become increasingly difficult to acquire, 
p and q assume values close to one and it can be proven that 
in the limit: 
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a ( 1-p) on 



( 11 ) 



and 



13 (l-q)oa 



( 12 ) 



If the attrition rate coefficient a is set equal to a(1 - 
p) and the coefficient b is set equal to 0(1 - <3)» then the 
two equations above are identical to Lanchester's equations 
for area fire. 



likelihood function for a heterogeneous aodel, however since 
the DAMSTAC model assumes homogeneous forces, the necessary 
equations have been reduced to a simpler form than that 
derived by Clark. By differentiating the logarithm of the 
likelihood function with respect to p, the following equa- 
tion is derived: 



In his doctoral dissertation, Clark developed the 




where 



Cj. = 0 if casualty k is RED 
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1 if casualty ic is BLOB 



ffl = number of BLUE survivors prior to Icth casualty 
k«i 

n, = number of RED survivors prior to kth casualty 

t, “ t, , = time between casualty !c & casualty k-1 
k f k-l 

Likewise, by differentiating the logarithm of the likelihood 
function with respect to a » and by solving for a , a second 
equation in terms of ^ and S is produced: 

“’E(1 - - -V. ) 

Differentiating with respect to q and 0 produces equivalent 

A A 

type equations which will allow the estimates for g and 0 to 
be determined in the same manner as follows. 

Generally, the two equations (13) and (14) can now 
be solved simultaneously to determine the values for p 
and a . Unfortunately, this is easier said than done. The 
complexity of these equations (more so in the heterogeneous 
case) does not allow explicit equations for p and a . The 
remainder of this chapter will discuss- a method for deter- 
mining the estimates of p and oi from these equations and a 
way to extract the necessary data from a high resolution 
Monte Carlo simulation run to compute these estimates. 
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E. ESTIMATING PAHAMETEBS FROM A HIGH RESOLUTION MODEL 
As mentioned before, equations (13) and (14) must be 
solved simultaneously to produce values for a and p. Gener- 
ally, these equations are too complicated to produce such 
values easily so other techniques must be used. The method 
suggested by Clark is to enter a value for p in equation 
(14) and solve for a . This value and the selected value 
for ? are then put in equation (13) and the equation is 
evaluated. Notice that both values are estimates since we 
are computing the maximum likelihood estimates of p and a • 
The values of p and a which give a value for equation (13) 
closest to zero are the best maximum likelihood estimates 
for p and a . If equation (13) equalled zero exactly, the 
true values of the maximum likelihood estimates would be 
known. 

In using these equations, the input received from the 
high resolution, Monte Carlo simulation must include: 
the time between casualties; 

the number of RED and BLUE survivors prior to the kth 
casualty ; 

the total number of casualties for RED and BLUE for the 
battle ; 

the type of casualty (RED or BLUE) of the kth casualty. 
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In Clark's theory for a homogeneous force battle, the 
conditional attrition rates and the probabilities that a 
target is not acquired are assumed to be constant over the 
entire battle. This means that equations (13) and (14) are 
summed over all casualties in the battle and the estimates 
for p and a , which are determined from them, are also 
assumed to be constant over the entire battle. 

In applying this theory to STAR output, the authors 
obtained questionable results for the values of the parame- 
ter estimates. The estimates, p and q, were computed to 
have values of 0.99 and 0.93 respectively. These values 
imply that both BLUE and RED forces experienced difficulty 
in acquiring targets. Analysis of the Fulda Gap terrain and 
consideration of the respective locations of the opposing 
forces support the estimated values for p and q. However, 
the stated value for p, 0.99, is a truncation of the actual 
value computed. If proper rounding procedures had been 
used, p would have been set equal to 1.0, which means that 
the RED forces would not detect any BLUE targets. Nondetec- 
tion implies that no BLUE attrition will occur. But this 
contradicts the STAR battle which did, in fact, produce BLUE 
casualties amounting to about two-thirds of the total BLUE 
force. 
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The estimates, aand j3, were computed to be 2.4 and 0.03, 
respectively. The relatively large value of a implies that 
when the RED forces are able to detect BLUE targets, the 
BLUE force will suffer heavy casualties. There appears to 
be an unusually large difference between the values for a 
and ^ which are not reflected by the values for p and g. 

This disparity is most likely attributable to the restric- 
tion that all paramete;c values be held constant over the 
entire battle. Experimentation with different parameter 
values in DAHSTAC indicates that a value of about 0.40 more 
closely approximates the attrition observed in the STAR bat- 
tle. Experience with these computations reveals that this 
is not necessarily the case and that a more likely assump- 
tion is that the parameters are constant over different time 
intervals of the battle. In general, a simple battle can be 
divided into phases where the intensity of the battle is 
constant within each phase but different among the phases. 
Basically, the phases are: 

movement to contact (low intensity) 
initial contact (mid intensity) 
main battle (high intensity) 

withdrawal from contact or annihilation (mid/low inten- 
sity) 
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Moving from one phase of the battle to the next should cause 
a change in conditional casualty rates ( aand/? ) as well as 
the probability of not acquiring a target (p and g) . 

The question now becomes how to determine the breakdown 
of the battle by phases using the output from a high resolu- 
tion, Monte Carlo simulation. That is, where do the changes 
in phases occur in the simulation? One possible method is 
to look at the number of casualties as a function of time. 
Another is to look at the number of shots fired as a func- 
tion of time (see Figure 11) . In both cases, time is seg- 
mented into time intervals (eg. every 100 seconds) and a 
histogram of the casualties or shots is plotted for each of 
these intervals. Having done this for several battles, the 
results identify fairly clear choices for the end of one 
phase of the battle and the beginning of the next phase. 
However, more analysis of these two methods reveals that 
they might not truly identify the phases of the battle and 
that an alternate method should be considered, primarily one 
based on range. This concept is discussed in detail in 
Chapter VI. 

Using the number of casualties per time interval, a dis- 
rortei picture of the intensity of the battle might be 
obtained. The interaction of the probabiliry of acquiring a 
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Figure 11: Histogram of Shots per Time Interval 

targez with the attrition rate of a firer aighr offset each 
other and produce results which conflict with the true bat- 
tle intensity. In other words, if the acquisition probabil- 
ity is high but the attrition rate is low, the number of 
casualties in a time interval mighr indicate a low inrensity 
phase of the battle when, in fact, it is a high intensity 
phase. Therefore, this method should be used with caution, 
if at all. The method of using the number of shots per time 
interval has similar faulrs. Twenty shots by one firer in a 
rime interval is much different in terms of intensity than 
one shot by twenty firers. Again, this method should be 
used with caution. 
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The question is then: what method should be used? One 
method, which is used by the Combined Arms Combat Develop- 
ment Activity (CACDA/CASAA) , is measuring the number of 
shots per firer per time interval. ^ Although it also has 
some faults, it does provide a better indication of the 
intensity of the battle. It is, for the most part, indepen- 
dent of the conditional casualty rates, although the number 
of shots would indicate a dependency on target acquisition. 
This dependency is reduced by averaging the number of shots 
over the number of firers. In effect, this measure can be 
considered to be unbiased by the changes in conditional 
casualty rates and probabilities of not acquiring targets. 
Thus one is able to get a feel for the intensity of the bat- 
tle from the histogram of the number of shots per firer per 
time interval. This is not necessarily scientific reason- 
ing, but it seems to work. 

Once the phases of the battle are delineated for a high 
resolution, Monte Carlo simulation run, the next step is to 
determine how to treat the *dead' time from the beginning of 
a phase to the first casualty and from the time of the last 
casualty to the end of the phase. For phases in the middle 

^Conversation with Mr. Ernest Boehner from CASAA on 5 
December 1981 
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of the battle, each of these two times is a portion of the 
time between two casualties and should thus have a bearing 
on our maximum likelihood estimates from equations (1 3) and 
(14) . 

To resolve this situation, the time of the last casualty 
within a phase is used as the end of the phase. For exam- 
ple, if a phase is supposed to end at 1200 seconds and casu- 
alties occur at 1140 and 1225 seconds, the phase will be 
terminated at 1140 and the new phase begun at 1141, This 
also allows this time between casualties to be accounted for 
in the next phase. This method produces some variation in 
phase lengths when running several replications of the bat- 
tle, but this variation is minimal, especially since the 
beginning and end of phases are being judged from something 
less than a scientific method. 

Several replications, as many as affordable, should be 
used to compute the maximum likelihood estimates for the 
four parameters. By using a large sample population (the 
casualties) , a more accurate estimate of the parameter can 
be obtained. Likewise, when the heterogeneous model is 
used, a larger sample population is desirable, since the 
number of engagements between two particular weapons systems 
will be much less than the number of engagements in the 
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homogeneous case. This could produce values that are 
statistically insignificant. 

A postprocessor to the high resolution model was devel 
oped to make the numerous calculations required to produce 
these estimators. The postprocessor uses an optimization 
routine to determine the maximum likelihood estimates from 
the given equations. Detailed information on the postpro- 
cessor is provided in Appendix E. 
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17. ALTERNATE F UNCT I ONAL FORM S 



Different military situations have been hypothesized to 
yield different functional forms for Lanchester-type equa- 
tions. The term "Lanchester-type equations" refers to any 
differential-equation model of combat and the term "func- 
tional form" relates to the structure of the attrition rate, 
i.e., how the attrition rate equation is written and which 
parameters are considered explicitly in the equation. 
Henceforth, a standard shorthand notation will be adopted to 
refer to the functional forms used in homogeneous-force Lan- 
chester-type combat models. For example, Lanchester's equa- 
tions of modern warfare are expressed as F|F functional 
forms (or, simply PjF attrition) because they are propor- 
tional to only the number of enemy firers. Taylor [Ref. 7] 
provides a good explanation of this convenient shorthand 
notation. 

Generally, the classical Lanchester functional forms for 
expressing the attrition rate for forces using aimed fire or 
area fire are too simple to be used "as is" in a detailed 
analytical combat model. The assumptions required for these 
basic equations to hold are too limiting and do not allow 
for all of the processes that are known to occur in combat 
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to be considered in the force attrition model. Some of the 



more common functional forms for attrition rates that have 
been considered and used in Lanchester-comba t models are 
shown in Table 2. 

Of special interest is the T|T functional form. This 
functional form is particularly useful in describing attri- 
tion in target rich environments. Peterson [Ref. 8] hypoth- 
esized that the equations, 

dx/dt = -ax and dy/dt = -by, 

characterize the early stages of a small unit engagement in 
which the vulnerability of a force dominates its ability to 
acquire enemy targets. It is recommended that Peterson’s 
T|T functional form be considered in the heterogeneous model 
described in Chapter 71 to model attrition in the initial 
phase of the battle. The remaining sections of this chapter 
will discuss several not-so-common functional forms that 
have been formulated and used in ’'recently” developed combat 
models . 

A. HELHBOLD-TYPE EQUATIONS 

An alternative functional form for homogeneous-force 
attrition rates was proposed by R. Helmbold. Helmbold 
believed that when opposing force sizes were grossly une- 
qual, the larger force experienced "inefficiencies of scale" 
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Table 2: FaNCTIONAL FORM NOTATION 



FUNCTIONAL FORM DIFFERENTIAL EQUATION 



F| F 


dx/dt 




-ay 


dy/dt 




-bx 


FT| FT 


dx/dt 




-axy 


dy/dt 




-bxy 


F|PT 


dx/dt 


s 


-ay 




dy/dt 


= 


-bxy 


TIT 


dx/dt 


s: 


-ax 


dy/dt 




-by 


(F+T) 1 (F+T) 


dx/dt 


= 


-ay - cx 


dy/dt 




-bx + dy 



in producing casualties. Limitations on space available for 
maneuver^ reduced reaction time, increased command, control, 
and communication problems, etc. result in the reduction of 
casualty producing effectiveness. These limitations may 
well prevent the larger force from using its full destruc- 
tive capability against the smaller force. There is no 
definition for ’’grossly" unequal force sizes, but the fol- 
lowing hypothesis is proposed. The validity of using Helm- 
bold’s equations in a combat model is directly related not 
only to the force ratios (i.e., x/y or y/x) but to the rela- 
tive size of the twc opposing forces with respect to their 
surrounding terrain. For example, a 5 to 1 force ratio may 
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well be a situation where Helmbold's equations apply if both 
of the opposing force sizes are relatively large with res- 
pect to the terrain available to them for movement, conceal- 
ment, etc. Conversely, a force with a 30 to 1 force ratio 
may not exhibit ineff iciences of scale if the opposing force 
sizes are small in relative sizes and with respect to their 
surrounding terrain. Consequently the decision to use Helm- 
bold's functional form for attrition rates rather than some 
other form (like Lanchester's equation for modern combat) 
really depends on the tactical "situation”. 

Helmbold's functional form introduces a modification 
that accounts for the reduced fire effectiveness of the lar- 
ger force. In mathematical form his equations would read: 
dx/dt = -a g (x/y) y with x(0) = x^ 

dy/dt = -b h(y/x) x with y(0) = Yq 

where a and b represent time-dependent atrrition-rat e coef- 
ficients, and g (x/y) and h (y/x) represent functions that 
modify fire effectiveness of individual x and y firers 
respectively. According to Helmbold's hypothesis the two 
functions g (x/y) and h (y/x) must satisfy: 
g(1) = h(1) = 1 
g(u) = h(u) = E(u) 

B(u) is a strictly increasing function of u. 
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So, the generalized Helmbold-type combat equations can be 
written as: 

dx/dt = -a E (x/y) y 
dy/dt = -b E(y/x) x 

For the case in which E(u) = u with c a constant 
greater than or equal to zero, the following equations for 
Helmbold-type combat are obtained: 



dx/dt = -a 


t« w 

(x/y) y 


with 


x(0) 


dy/dt = -b 


(y/x)‘"'^ X 


with 


y(0) 



H is commonly called the "Weiss parameter" with W = 1-c. 

The above equations represent the special exponential 
case of Helmbold combat. If a(t) and b(t) represent constant 
attrition rate coefficients denoted by a and b respectively, 
then, when W = 1, the above equations reduce to Lanchester’s 
Equations for Modern Combat. When W = 1/2 Lanchester’s 
Area-fire Equations are obtained. 

B. HOW TO CHOOSE THE RIGHT FUNCTIONAL FORM 

A very logical question that readers of this study may 
ask at this time is, "Now that I've seen some of the Lan- 
chester-type attrition rate functional forms, which one 
should I choose for my analytical model?" The answer to 
this question is that there is no answer to this question. 
The choice of functional form is purely arbitrary and is 
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usually left to the discretion of the model builder. The 
following steps, however, may help in this choice. 

1. Carefully consider the scenario that the combat model 
is intended to portray. The scenario may well indicate 
that some particular subset of the functional forms is 
preferable. 

2. List the assumptions that must be made for each of the 
functional forms in the subset to hold and see which 
assumptions are acceptable to you and your model and which 
are not. 

3. Observe the parameters involved in each functional 
form and eliminate those forms for which the parameter 
values or estimates are unobtainable or for which you are 
not willing to plug in arbitrary values. 

At this point, the final selection of functional form is 
purely one of subjective judgement. 

There is no requirement, however, that a model builder 
must accept one of the already developed functional forms. 
There is no doubt that other forms can be formulated, but 
the formulation will require extensive research, analysis 
and work. A good check to see if an attrition rate equation 
is reasonable is to insure that the Lanchester classical 
equations are special cases of the new functional form. 
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Once the final functional form has been selected or 



formulated, it may be enhanced as necessary to reflect more 
of the complexity and realism of combat. Some additional 
operational factors that can be considered are: 

Onit breakpoints and battle termination criteria; 

Suppression ; 

C3; 

Reinforcements or replacements. 

The main assumption that allows these operational factors to 
be "added” to the attrition rate equation is that synergis- 
tic effects are not considered in Lanchester-type force-on- 
force combat models. No synergistic effects imply that the 
Principle of Superposition holds, therefore, the effects of 
the operational factors are independent and additive. 
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V . FUTURE ENH ANCEMEN TS 

There are numerous enhancements to the analytical model 
that can be made. The enhancements that are recommended, 
for the most part, are not aimed at improving the inner 
workings of the current model, but at building successive 
models that are more detailed, require fewer limiting 
assumptions, and are generally more complicated. 

For follow-on models that retain the homogeneous force 
concept, it is recommended that the model have the ability 
to allow the user to choose alternative functional forms for 
the force attrition rates. This should be done to better 
describe the situation of the battle at a particular time. 

A likely method would be to change from one functional form 
to another during rhe battle in order to more realistically 
portray the attrition occurring then. This enhancement 
would require considerable effort in developing the maximum 
likelihood estimates for the various parameters used in each 
alternative functional form, additional work with the STAR 
postprocessor, and significant change to the current DAMSTAC 
model. Questions which arise from this suggestion are: 

Which forra(s) should be chosen? 
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How do you know when to change functional form? 

These questions are partially answered in Section I7-C. 
Further research into this area is required. 

The movement subroutine of the model needs to be modi- 
fied to more realistically portray the movement of units on 
the battlefield. Even though a relatively low-resolution 
time-step operation is occurring, the resolution of the 
movement of units can be refined to vary the rate of move- 
ment from one time step to the next, if necessary. Cur- 
rently, a change in rate of movement occurs only at the 
initial contact between forces and at breakpoints. The rate 
of advance of the RED forces should be allowed to vary based 
on decision logic which considers attrition rate, terrain, 
etc. BLUE forces should be given the capability to move to 
alternate or fall back positions when the batxle situation 
dictates. Decision logic must be enhanced to allow BLUE 
forces to move to alternate or fall-back positions. It 
should be noted that this "pace of battle" question has 
never been fully explained or understood, however, the move- 
ment routine logic in DAMSTAC should parallel that of STAR. 

The concept of units having different formations is 
ignored in the current model. Units are identified by a 
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point location only. some work needs to be done to allow 
the units to take on different formations for road marches, 
attacks, etc. In order to do this, additional attributes 
such as frontage, depth, orientation and geometric configu- 
ration should be considered to describe the formation of the 
unit. a postprocessor for STAR could take the individual 
tank locations from the STAR data base and determine the 
above mentioned attributes for input into the DAMSTAC model. 
Attrition rates on both sides must be adjusted for the vari- 
ous formations of an attacking unit. 

The inclusion of indirect fire, suppression, smoke, 
intelligence, and close air support are all important to 
making the model a usable analytical tool. Additions of 
communications and logistics, primarily ammunition and 
petroleum, oil, and lubricants (POL) , are also important 
enhancements that must be included in the model. 

As was seen in Section III-E-3, P^^,the probability that 
a target is visible and acquired in the steady state is 




Both of the factors represent a particular element of the 
attrition process, one the terrain affects on line of sight. 
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and rhe other# the acquisition process. It should be appa- 
rent that these two factors can he separated from each other 
with the only common variable, p the rate of transfer from 
the visible state to the invisible state. With this factor- 

V 

ization, if the line of sight parameters for a particular 
terrain type are know.n and the attrition rate coefficients 
have been determined, the target acquisition capability of a 
weapon system against another weapon system can be deter- 
mined. This would allow the target acquisition parameter to 
be used for other terrain types to compute attrition rate 
coefficients without having to estimate the coefficients 
from a high resolution simulation. 

The concept works in this manner. Suns are made on the 
high resolution model (STAR) using a particulcir terrain. 

The attrition coefficients are determined using the maximum 
likelihood estimates described in Section III-D. with pre- 
determined data for the terrain coefficients, the line of 
sight parameter can be factored out of the attrition rate 
coefficient, A. This would leave the value of the acquisi- 
tion factor which should theoretically be constant for the 
firer-target type combination used in the scenario. With 
this value and the terrain coefficient for other types of 
terrain, the attrition rate coefficient can be determined 
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for a new scenario with the same firer-target combination, 
without running the high resolution model again. This 
method assumes that the terrain affects and the target 
acquisition process are independent of each other. The 
impact of this concept is enormous in terms of time and cost 
saved from not having to make multiple runs of the high 
resolution model. 

One concept which is not strictly an enhancement but is 
an important consideration is the idea of using simulated 
combat data from a source other than STAR. Specifically, 
data generated at the National Training Center from actual 
field exercises could conceivably be used to compute attri- 
tion coefficients in much the same manner as STAR. What 
this provides is not only an alternate source of data, but 
also a different model, since it would include many human 
factors that would not likely be accounted for in STAR. 
Depending on the intended use of the model, one source or 
the other may be desirable. 

It is necessary that a progression from homogeneous 
forces to heterogeneous forces be made if DAMSTAC is ever to 
be really useful. This enhancement entails a significant 
amount of research, model building, and statistical analysis 
in order to implement the concept of heterogeneous forces 
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into the analytical model. Target priorities, fire 
allocation algorithms, and a library of attrition rate coef 
ficients must be developed for various scenarios and force 
mixes. Changing the model to a heterogeneous type requires 
some changes in the equations (13) and (14), and thus some 
significant modifications to the postprocessor for the 
model. 
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VI . H ETER OGENEOUS MQD^ 



DAMSTac, in its current configuration, models combat 
between two opposing homogeneous forces, but actual combat 
consists of many different weapon-system types operating 
together as combined arms teams. For example, there may be 
infantry, tanks, antitank weapons, artillery, etc. on each 
side. A necessary extension or modification, therefore, to 
DAMSTAC is the consideration of combat between heterogeneous 
forces. This chapter briefly indicates how the basic ideas 
and methodologies used previously to model combat attrition 
for homogeneous forces can be extended and adapted to the 
heterogeneous force structure. 

A key assumption that must be made in order to model 
combat between heterogeneous forces, where Lanchester-type 
differential equations are used, is that the attrition 
effects of various different enemy weapon systems types 
against a particular friendly target type are additive. In 
other words, synergistic effects between weapon system types 
are not considered [Ref. 9]. Also inherent in the modelling 
of heterogeneous force combat is the need to introduce an 
allocation factor. This allocation factor may be defined as 
the fraction of the firers of a particular type that engage 
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a particular target type. The allocation factor concept 
will be discussed in greater detail later in this chapter. 

One theoretical concept for modelling the attrition pro- 
cess between two opposing forces is to consider attrition as 
a nonstationary Markov process where the states of the pro- 
cess are defined as the numbers and types of surviving com- 
batants and the ranges between opposing combatants. The 
transition probabilities between the Markov states, under 
this concept, depend on the ratio at which surviving combat- 
ants or weapon systems inflict casualties on a particular 
type of opposing combatant. Thus the nonstationary Markov 
process model provides a means of relating the sequence of 
states and the times of the state transitions observed dur- 
ing a battle to attrition rate coefficients. In turn, this 
allows for the estimation of attrition rate coefficients on 
the basis of the observed sequence of states and transition 
times. Andrighetti [3ef. 10] developed such a ncnstationary 

Markov model that related weapon system effectiveness to the 
time sequence of casualties observed in a two-sided, hetero- 
geneous force land combat simulation. There is a flaw, how- 
ever, in Andrighetti' s work that readers should note. He 
implies that his methodology is easily extended to time-de- 
pendent attrition rate coefficients. This is not true. 
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Andrighetti *s thesis does contain a good discussion on 
the theoretical background which allows for the maximum 
likelihood estimation of parameters used in attrition rate 
equations. An important observation made by Andrighetti in 
his thesis is that the difference-differential equations 
obtained from the Chapman-Kolmogorov equations of the 
nonstationary Markov model are susceptible to a general 
solution only in special cases as Gordon Clark showed in his 
doctoral dissertation. Clark used the Chapman-Kolmogorov 
equations to develop expressions for time state probabili- 
ties and the expected survivors in a heterogeneous force 
battle. 

The extended version of the DAMSTAC model which shall be 
called DAMSTAC II should be a heterogeneous force analytical 
model that uses STAR output to compute parameter values. 
These parameter values should be allowed to change as the 
distance or range between opposing forces vary during the 
battle. In DAMSTAC, an attempt was made to allow the param- 
eter values to change with respect to time. However, it was 
discovered that the maximum liklihood estimators developed 
for the critical parameters in the attrition rate equations 
could not be used for time dependent variables. Two feasi- 
ble alternatives that allowed for variable parameter values 
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were considered. Both alternatives are recommended for 
future evaluation. The first alternative looks at dividing 
the battle into phases based on battle intensity as measured 
by the number of shots divided by the number of live firers. 
(Refer to Section III-E) The second alternative for divid- 
ing the battle into phases is based on the range between 
combatants. Here also the parameters can change from phase 
to phase. This technique will work assuming that the com- 
batants on each side move as a single force and that units 
are located by their centers of mass. For example, a battle 
can be divided into three distinct time intervals or phases 
of combat: an artillery prep phase during which all combat 
vehicles and weapon systems, other than artillery, are not 
participating; a moderate to long range battle phase; and, a 
terminal shorter range battle phase. Parameter values can 
then be computed, using procedures to be discussed later, 
for each battle phase. This second alternative is used in 
the COMANEW model and is the one recommended for DAMSTAC II. 

A. THE ATTRITION-RATE FaNCTIONAL FORM 

The attrition rate functional form that will ultimately 
be chosen for the heterogeneous force model is, of course, 
the prerogative of the modeller. The authors of this the- 
sis, however, would like to recommend two related attrition 
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rate functional forms for consideration. The first, which 



shall be labeled as Functional Form 1 (FF1) , is based on the 

assumption that during the course of a battle fire is 
directed at the highest priority target surviving and 
acquired. Also it is assumed that the relative priority of 
target weapon types, as considered by each firer, remains 
constant throughout the engagement and that this relative 
priority is identical for each firing weapon type. Func- 
tional Form 1, developed by Gordon Clark for use in the 
COMAN model, can be written as follows: 



type j during phase i of the battle given strengths of m 
and n. 

U) 

f^ (m,n) = analog of f^^. for RED weapons of type k. 




(16) 




where 

(m,n) = conditional casualty rate of BLUE weapons of 
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M 




= the number of surviving RED weapons of type k 



m- = the number of surviving BLUE weapons of type j 



= kill rate of a BLOE firer of type j for an acquired 
RED target of type k in the i'th time interval. 



“ijic" kill rate for a red-f irer/blue-target combination 
corresponding to /S;^ 



qj = the probability that a specific RED weapon is unac- 
quired by an individual BLUE firer in the i'th phase of 
the battle. 



p; = the probability for a specified bl ue-target/red-firer 
combination corresponding to g; . 



mj = the number of surviving BLUE weapons of higher prior- 
ity than a type j weapon. 
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= the number of surviving RED weapons of higher prior- 
ity than a type Ic weapon. 

The reader should note that since the assumption is made 

that a firer will engage the highest priority target he can 

acquire, the probability of a BLUE firer engaging a RED 

weapon of type Ic is qj'^(1- Similarly, the probability of 

/ 

a RED firer engaging a BLUE weapon of type j is . 

These two probabilities may be regarded as the allocation 
factors for Functional Form 1. 

All of the critical parameters necessary to compute the 
conditional casualty rates, f*{^j (m,n) and f"J^ (m,n) , can be 
obtained either directly or indirectly (using maximum like- 
lihood procedures) from the output of a high resolution 
simulation such as STAR. 

If one is not willing to assume that target priorities 

remain constant throughout the battle then FFI will not be 

sufficient. The terms in Clark’s functional form which 

result from the assumption of constant target priorities are 

qj'’ and p-^i where qj'^is the probability that all surviving RED 

weapons of higher priority than type k have not been 

/ 

acquired and p^^ is the corresponding probability for type j 
blue weapons. 
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One way to extend Clark's formulation in order that the 
fixed priority may be dynamically overridden is by replacing 
the terms q. with0j^ where is defined as the probability 
that BLUE weapon type j engages RED weapon type k, given 
that BLUE weapon type j acquires at least one RED weapon 
type k. Note that 

^j^(1-q*'*^) = P(j engages k|at least one k acquisition) 

P (acquire at least one k) = P(j engages k) 



Functional Form 2 (FF2) may be written as; 



K\ = E (1-pP) 



m: 



(17) 



(m,n) ) 



(18) 



where is the RED firer BLUE target probability corre- 
sponding to 9J^^ and all other parameters are as defined for 
FF1. Just as p,q,Q!,j3 are assumed constant over the entire 
phase of the battle, 5 and 9 , for each j and k, are also 
assumed constant over the entire interval. 
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Serious consideration should be given to incorporating 
FF2 or some similar functional form into DAMSTAC II because 
STAR does allow target priorities to change during a battle. 
Target priorities in STAS are based on the types of targets 
on the firer*s detected list and the type of ammunition the 
firer has available for use. 

It may also be useful to use a different functional form 
during the initial phase of combat where attrition rates are 
more a function of the exposure or vulnerability of targets 
than they are a function of the firer *s capability to 
acquire targets. Clark recommended the use of Peterson's 
Logarithmic Law to model attrition during the initial phase 
of the battle. [Ref. 11] 

B. MAXIMOH LIKELIHOOD ESTIMATORS 

The ability of the DAMSTAC II model to provide insight 
into the interactions represented in the STAR simulation is 
based on the estimation of conditional casualty rates ( a 
and )3 ) and other parameters from STAR data. These estimated 
parameters should be constant within a phase of the battle, 
but independent and unrelated to results in neighboring time 
intervals. The independence of parameter values between 
time intervals implies that the estimators for a specific 
time interval, or phase of the battle, can be defined by 



85 



analyzing a sample of data from the high resolution simula- 
tion observed during that single -time interval. 



The maximum likelihood estimators (ULS's) for the 
unknown parameters p and g used in FF1 were developed 

by Gordon Clark. It may be of interest to the reader to 
note that simplified versions of these ML2*s were used in 
the DAhSTAC model for homogeneous forces, however, since we 
are now concerned with heterogeneous forces, the ’’una- 
bridged" versions are presented. The maximum likelihood 
estimate for or as a function of p is: 



The estimators obtained from the above equations result 
from the data from a single replication of the combat simu- 
lation. done precise estimates are obtainable by using a 
larger sample consisting of a number of observations during 







for > 0 



( 19 ) 



0 



Otherwise 



where the estimator for the parameter o; for a given p 

The MLE for p is found by solving: 



9ln L(d) 
3P 





(20) 



) . 



(niy^ -(myj )) = 0 

The aquations for (q) and q are similar. 



>fcl jml \9i 
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the interval in question from a number of independent 
replications of the STAR simulation. This modification 
entails an additional summation , where x = the total 

X 

number of replications) in each of the estimating equations. 
The data needed for the computations of the MLE's include; 
y = total number of casualties 

ty = elapsed time since the beginning of the interval 
until the y*th casualty 

ty^,= the end of the battle phase or the elapsed time since 
the beginning of the phase until the battle ended if it is 
less. 

m^ = the number of BLUE survivors of type j just prior to 
the y'th casualty 

n^,^ = the number of RED survivors of type Jc in the battle 
jusr prior to the y'th casualty 

C>f = 1 if the y'th casualty is BLUE ; 0 otherwise 

t>f = firer weapon type for the y'th casualty 

= target weapon type for the y'th casualty 

myj = the number of BLUE survivors of higher priority than 

-YP® j just prior to the y'th casualty 
^ b 

®'ri ~ S 

ttk I 

n^ = RED analog of niYi 
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the number of BLOB weapon-type j casualties due to 




RED weapons of type k 

One technique of computing the MLE*s for a and p is to 
numerically determine p using equation (20) and then substi- 
tute this value into equation (19) for j = 1,2,...,b and for 
k = (b = the total number of BLUE weapon types and 

r = the total number of RED weapon types). A recommended 
technique for numerically solving for p, the estimate for p, 
can be found in Clark *s work on the Dyntacs combat model 
[Ref. 12]. This technique is similar to the one developed 
for the homogeneous case discussed in Chapter III. Basi- 
cally the procedure is as follows. Since p represents a 
probability, it must satisfy the inequality 0<p<1. There- 
fore, the function ^ Ln L(p) can be evaluated at points in 
the interval 0<p<1 until two points ,p, and p^ , are found 
such that the product of Ln L (p, ) and Ln L (p^) is less 
than or equal to zero. Since Ln L(p) is a continuous 
function,-^ Ln L(p) =0 for some p satisfying p, ^ p ^ Pj . 
This interval is reduced by evaluating Ln L(p) at the 

midpoint of the interval (p, ,Pj). Repeating this process 
converges to a solution for Ln L(p) and this solution is 
a local maximum. This technique is similar to the bisection 
method for solving nonlinear problems. 
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with respect to Functional Form 2 , the introduction of 
the 9jj^ allocation terms adds considerable complexity to the 
maximization of the likelihood function. One way to circum- 
vent this problem is to compute a "best guess" for the *s 
from STAB samples (for the given interval) . This best guess 
can thenbe substituted into the ©j,^ *s prior to maximization 
of the likelihood function, thereby reducing the©jj^'s to 
constants. The maximization can then be performed, as in 

I 

Clark's original equations, with respect to p, q,Q; , /3 . 

C. THE MODEL 

The computational procedures that will constitute 
DAMSTAC II must be developed by the modeller who chooses to 
carry on the work initiated in this thesis. As a starting 
point, however, the techniques used by Clark in his COMAN 
model would be excellent choices. The adequacy of Clark's 
procedures in DAMSTAC II given a particular scenario and 
force mix, to predict attrition in a STAR simulation with 
the same scenario and force mix is questionable and would 
have to be tested. Detailed comparisons and analysis of the 
output from both models will be required to answer the ques- 
tion of model comparability and modifications to the com- 
puter code and/or computational procedures in DAMSTAC II 
must be accomplished as necessary. The computational cycle 
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employed in COHAN is based on the assumption that the dis- 
tribution times between casualties in the high resolution 
simulation is exponentially distributed. It uses Honte 
Carlo procedures to generate the time until the next casu- 
alty and also to find the casualty weapon type. 

The following steps constitute the computational cycle 
recommended for DAHSTAC II. 

Step 1. The first step in the cycle sets rhe time equal 
to zero and sets the battle phase, designated by the letter 
i, equal to 1. 



Step 2. Next the total conditional casualty rate, 

\‘(m,n) is computed from the equation 
\'Cm,n) =1^^ f“J.(m,n) + for i= 1,2,---^ I (21) 

where f^ (m,n) and f^ (m,n) are the conditional casualty 
rates in phase i of the battle, given force strength m and 
n, for BLUE weapon type j and RED weapon type k, respec- 
tively. \* (m,n) is the parameter of the exponential distri- 

bution describing the times to the next casualty. 

Step 3. Now, using the exponential distribution with 
parameter X*(m,n), one Honte Carlo's for the time interval 
to the next casualty. Designate this time interval as t'. 

If the next casualty occurs in another time interval, ie. if 
t + t' > t; where t*, represents the end of the i'th battle 
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phase, then i is increased by 1 , the time is set equal to t 
and the program is directed back to step 2 where \''*‘(m,n) is 
selected as the distribution parameter. If, however, t ♦ t 
< t; then the next step is to Monte Carlo for the casualty 
weapon type. 

Step 4. The weapon type for the next casualty is repre 
sented by the random variables and which are defined 
as 

Cjj = j if the casualty weapon is BLUE weapon type j 
0 if the casualty is RED 

Cp = k if the casualty weapon is RED weapon type k 
0 if the casualty is BLUE 

The Monte Carlo procedure uses the conditional casualty 
weapon probability distributions for and Cp . These 
probabilities are computed by 

P' (Cb = j|C = 1,m,n) = f^^ (m#n) 
p* (Cp — k|C— 1,m,n) ~ f (m , n ) / A^in^n) 

for i= 1,2,..., I and when variable C is set equal to 1 if 
a casualty is determined to occur in interval i from the 
previous steps. The probability that the casualty is of a 
particular type is computed for each RED and BLUE type 
weapon. The sum of these probabilities should equal 1. 
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step 5. The next step is to Monte Carlo using the Uni- 
form (0,1) distribution to determine the casualty type. An 
example may help in understanding the above procedure. Sup- 
pose there are two types of BLUE weapons and two types of 
RED weapons. Assume the conditional casualty probabilities 
for each weapon are computed using the equations shown pre- 
viously and they are 0.25 for each weapon type, b1, b2, r1, 
r2. If one orders the weapon types in some arbitrary fash- 
ion such as b1, r1, b2, r2, and then accumulates the prob- 
abilities on the interval (0,1), one can assign an interval 
to each weapon type. The example would yield the following 
int srvals: 

b1 - (0,0. 25) ; 
r1 - (0. 26,0.50) ; 
b2 - (0.51,0.75) ; 
r2 - (0. 76, 1 .00) . 

The uniform random number can then specify the casualty 
weapon type. If the random number is 0.60 then the casualty 
weapon type in our example is b2. 

Step 6. After the casualty weapon type has been found 
the next step is to decrement the number of that type weapon 
by 1. Now increase the value of t by t* , the rime to the 
next casualty. If t is greater than or equal to t^ which 
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represents the ending time of the last phase of the battle, 
or if either side is annihilated , the battle has ended and 
the program terminates. If t is less than tj then go to 
step 2. 

Observe that the above computational procedures are sto- 
chastic in nature and different replications of an engage- 
ment will likely yield similar but different results. If 
this method is chosen for DAMSTAC II, variance reducing 
replications for each engagement should be performed. 

D. THE STAR POSTPROCESSOR 

It is very important that anyone who attempts to develop 
the DAMSTAC II model understand and be able to manipulate 
the STAR simulation model. In the initial attempts to 
extract data from STAR, certain subroutines must be turned 
off. These routines include the SmoJce, Engineer, Close-air 
support, etc. whose effects on parameter estimates used in 
DAMSTAC II have not been thoroughly analyzed. The number of 
types of combatants played in the simulation will have to be 
regulared so that they do not overpower the capabilities of 
the analytical model. It is recommended that no more than 
three different types of weapon sysrems for each opposing 
force be attempted when initially formulating the heteroge- 
neous force model. 
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It is important that the finished version of the 
analytical model be as simple and transparent as possibl'e 
and not encumbered with parameter estimating algorithms and 
computations. This consideration plus the amount of data 
that must be extracted from the STAR output require the 
development of a STAR postprocessor program. As an illus- 
tration, if one considers a simple battle divided into three 
battle phases with only three types of weapon systems on 
each side, the DAMSTAC II model, as it has been described, 
requires the estimation of 114 parameter values and the com- 
putation of six attrition rates. 

The postprocessor at Appendix E, developed for the homo- 
geneous model, should provide a good foundation for the more 
complicated postprocessor needed for DAMSTAC II. 
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VII . CONCLUDING REMAR KS 

The concept of using differential equations to represent 
combat began with Lanchester’s classical equations and has 
developed into a very dynamic field of study. The theoreti- 
cal background of the fitted parameter method has been well 
documented by Clark and Taylor. However, the practical 
application of the theory lacks this same degree of discus- 
sion and documentation. This paper has focussed on some of 
the techniques of applying Clark's theory to an actual high 
resolution, Monte Carlo simulation. The problems are real 
and, as has been seen, the solutions presented are not 
doctrine. It is hoped that further investigation, elabora- 
tion and documentation of this modelling technique will 
occur in the near future. 

The choice of which functional form to use in a particu- 
lar analytical combat model is entirely up to the discretion 
of the model builder. In choosing an attrition rate func- 
tional form the model builder should consider the scenario 
portrayed in the model and the assumptions that are compati- 
ble with it. 

Finally, it is recommended that the entire program for 
the model be converted to the UCSD Pascal language on the 
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Apple II. The basic structure of Pascal is more conducive 
to combat modelling simulation than the structure of the 
BASIC language. The concept of entities and attributes 
found in SIMSCRIPT II. 5 can be closely emulated in Pascal 
while this is not as easily accomplished or clearly defined 
in BASIC. 

The modular structure of procedures and functions in 
Pascal is also an advantage in the enhancement process of 
model building. Modellers can use one procedure to repre- 
sent a particular action or phenomenon occurring on the bat- 
tlefield. Any action which is unable to be modelled can be 
left in the program as a dummy procedure for later work. 

Additionally, the Pascal language allows more descrip- 
tive variables to be used than the two-character variables 
used in Applesoft. This would make the program itself more 
transparent and easier to modify and update. 

The DAMSTAC model is clearly in its infancy stage. 
Through continued enhancement, this model can contribute to 
the ongoing work in the field of combat analysis. 
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APPM5IX A 



LIST OF VA RIAB LES 

A = casualty rate for RED attritting acquired BLUE targets 
AN = angle in radians measured counter-clockwise from due 
east made by the line between the current unit 
location and the next coordinating point. 

B = casualty rate for BLUE attritting acquired RED targets 
BE = the breakpoint for the BLUE force 
BM = movement criteria for BLUE units as a fraction of 
unit strength 

BF(I) = number of RED targets being fired at by BLUE unit I 
BP (I) = most recently passed coordinating point for BLUE 
unit I 

BR(I) = rate of movement for BLUE unit I in M/SEC 
BT(I,J) = attrition (casualties) caused by RED unit J on 
BLUE unit I in one on one battle 
BX(I,J) = x-coord of coordinating point J for BLUE unit I 
BY(I,J) = y-coord of coordinating point J for BLUE unit I 
BO (I) = current total BLUE force strength for time period I 

BO (0) = initial BLUE force strength 

B1 = number of coordinating points for BLUE (including 
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initial location) 

B1 (I) = initial strength of BLUE unit I 

DB (I) = number of casualties or amount of attrition for 
BLUE unit I 

DH = distance travelled by a unit in time increment DT, at 
speed RT 

DI = distance between current unit location and next turn 
point 

DR ( J) = number of casualties or amount of attrition for RED 
unit J 

DT = time increment 

DX = X component of distance between current unit location 
and next coordinating point 

DI = y component of distance between current unit location 
and next coordinating point 

DI = X component of distance moved in time increment DT, at 
speed RE (I) or BE (I) 

D2 = y component of distance moved in time increment DT, at 
speed SRtD or BE (I) 

a = number of BLUE units 

MT = max time of simulation 

H = number of RED units 

N$ = file name (used recursively) 
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p = probability that a BLUE target is nor acquired by 
a RED unit 
PI = 3. 14159265U 

PR = status of printer; 0 = printer is off 

1 = printer is on 

Q = probability that a RED target is not acquired by 
a BLUE unit 

RE = the breakpoint for the RED force 

RF (J) = number of BLUE targets being fired at by RED unit J 
RG = range between a BLUE unit and a RED unit 
RM = movement criteria for RED units as a fraction of 
unit strength 

RP(I) = most recently passed coordinating point for RED 
unit I 

RR(iy = rate of movement for RED unit I in M/SEC 
RT(J,I) = attrition (casualties) caused by BLUE unit I on 
RED unit J in one on one battle 
RX(I,J) = z-coord of coordinating point J for RED unit I 

RY(I,J) = y-coord of coordinating point J for RED unit I 

R0(l) = current total RED force strength for time period I 

RO (0) = initial total RED force strength 

R1 = number of coordinating points for RED (including 
initial location) 
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R1(I) = initial strength of RED unit I 
SB (I) = size of BLUE unit I 



SR (I) = size of RED unit I 



T = simulation 


time 








XB(I) 


= current 


x-coord 


for 


BLUE 


unit I 


XR (I) 


= current 


x-coord 


for 


RED 


unit I 


YB(I) 


= current 


y-coord 


for 


BLUE 


unit I 


YR(I) 


= current 


y-coord 


for 


RED 


unit I 
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APPEND I.X B 



NOTES ^ THE COMPUTER MO DEL 

1. Data which have been hardwired into the model are: 

A = 2.4 

B = 0.03 
BE = 0.3 
BT = 0 
DT = 30 
MT = 1500 
P = 0.99 

PI = 3. 141592654 
Q = 0.93 
RE = 0.3 

2. Once contact is made (range < 2700) , the rate of advance 
decreases to 2 meters/second. 

3. The option exists within the program to allow the move- 
ment of units based on a criteria of fractional force 
strength. When unit strength goes to 50% (movement crite- 
ria) , the unit will move to an alternate location (for BLUE 
units) or retreat along its path of advance (for RED units) . 
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A PPENDIX C 



utility PBOGRAMS 

This appendix provides a brief description of the util- 
ity programs associated with the Aggregated Model of STAR on 
the Apple II Computer. These programs provide the capabil- 
ity to the user to transfer the necessary data from the STAR 
model to this model with the minimum effort. 

The utility programs available are: 

Force Maker 
Results Reader 

Both programs are interactive with the user, who answers 
questions provided by the computer. These programs use the 
commands for the Apple II Disk Operating System (DOS) which 
allow them to read/write Sequential Text Files and Random 
Access Files onto the disk. Familiarity with these types of 
files and the DOS would improve the understanding of these 
utility programs. [Ref. 13] 

In the following program descriptions, the variable name 
used in the program is listed in parentheses after its 
description. 
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A. FORCE MAKER PROGRAM 



This program will write data files containing informa- 
tion about the two opposing forces. 

Force Maker produces a Sequential Text File containing: 
Force size (number of units) (M) 

The strength of each unit (SB (I) ) 

The number of coordinating points for movement of each 
unit (the first coordinating point is always the initial 
location of the unit) (B1(I)) 

The grid coordinates of each coordinating point for each 
unit (BX(I,J), BY(I,J)) 

The user has the option to make a file for the BLUE 
force or the RED force. After all the data has been 
entered, a data check routine will print the information for 
all units, so that a check for errors can be made. After 
the user has checked the data and identified the name for 
the file, the program will record the file on the diskette. 

A copy of the program is included as Listing B in Appen- 
dix F and a sample session using this program is attached as 
Figure 12, 
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3FI.UN FCRCE i-IAKER 
r':£NU: 

'* * ■ 2 T 



■ r'iAKF 



FOR 5L'J£ FCRCS3 
2 ^-^AKE FILE FCR REE FCRCE3 

O! •• T T 



■WHICH? i 

WHAT IS THE FORCE SIZE CF SLUE 3 
FCR LJNIT 1 

hCW MANY CCGRDIMATIMQ POINTS (IMCLUDIM6 START LGCATICM) ? 
ENTER '3RID CGCRDS FDR CGCRDIMATINB POINT ^ 1 

v_prinQn rc sT/i-) 

\ w W . W W W W *T i J 

Y-CGQRD IS 1 004j25 
FCR L-NIT C 

HGW MANY CGGRDIMATING POINTS (INCLUDING START LGCATIQM/ ? 

ENTER GRID CGCRDS FOR CGCRD I MATING !=-GINT * 1 
:<-CGCRD IS 547SS 
r-CGQRD IS 98120 

-CR UNIT 3 

HCW MANY COORDINATING POINTS (INCLUDING START LCCAT 

ENTER GRID CGCRDS FOR COORDINATING POINT # 1 
X-CGGRD IS 37-148 
Y-CGGRD IS 



UNIT FT X-CGGRD r -COORD 

1 1 341. S3 100s33 

2 1 347SS '98120 

3 . 37443 ^7893 

ANY CHANGES? (Y/N) M 
NAME CF FILE*' THESIS 

i MAKS FIl_E FCR SLUE FCRC 
Z MAi-E FILE FOR 



;-ED FORCES 



'~!j 



Figure 12: Sample Session of Force Maker Program 
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B. RESULTS READER PROGRAM 



This program reads the data file of time intervals and 
corresponding BLUE and RED force levels produced by the 
model. Results Reader produces two sequential text files: 
the first file with suffix ".COEFS” contains: 

Conditional Casualty rate, a (A) 

Conditional Casualty rate, /3 (B) 

Probability that a BLUE target is not acquired by a RED 
firer (P) 

Probability that a RED target is not acquired by a BLUE 
firer (Q) 

The second file with suffix ''.RESULTS" contains: 

The number of data points for BLUE force level versus time 
(N) 

The number of data points for RED force level versus time 
(N) 

Tima of the battle at time interval I (T) 

BLUE force level corresponding to time interval I (B0(I)) 
Tima of the battle at time interval I (T) 

RED force level corresponding to time interval I (R0(I)) 
The ".RESULTS" file contains some duplication of data for 
the following reasons. With the ".RESULTS" file written in 
this format, and free of additional data such as that in the 
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".COEFS” file, it can be used directly by the Apple Plot 
Program to produce a plot of the force levels versus time. 

A listing of the program is included as Listing C in Appen- 
dix ?. A sample plot of the DAMSTAC output data produced by 
the Apple Plot Program is included as Figure 13. A sample 
of tha output of DAdSTAC from Results Reader is included as 
Figure 14. 
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Figure 13: Sample Plot of Results Reader Program Output 
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APPENDIX D 

terrain model 

In addition to the basic computer model, work was done 
to input the STAR terrain data base into the Apple II compu- 
ter. This effort was done during the initial phase of the 
thesis work based on an erroneous assumption that the 
authors made concerning the functional form of the attrition 
rates. This assumption was that the terrain of the battle- 
field and line of sight computations must be explicitly 
played in the model. In the given functional forms for the 
attrition rates, line of sight and terrain are accounted for 
by the probability values of p and g. 

Before this assumption was found to be false, four util- 
ity programs, and a model that computed the elevations of 
the units were developed and are included in this document 
for possible future use. If further work is done with a 
model on the Apple II computer, the ability to input the 
STAR terrain into the Apple II may become an important part 
of the model. Examples of models that might use the STAR 
terrain in the Apple II are small aggregated models with 
functional forms that explicitly play line of sight, or 
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small, high-resolution models that are suitable on the Apple 
II computer. 

A. THE TERRAIN UTILITY PROGRAHS 

The utility programs that allow the STAR terrain to be 
input in the Apple II computer are listed below. 

Map Maker 
Hill Maker 
Woods Maker 
Forest Maker 

These four programs have two portions each: a file maker 
and a file reader. In general, the file maker portion 
allows the user to input a new file or update/change an old 
file. The file reader allows the user to print out on the 
screen or to a printer the selected file. 

For the Apple II combat model, a STAR terrain file (eg. 
FOLTERR file) is divided into its four main parts; 

the list of hill reference numbers in each grid square; 
the hill parameters; 

the list of forest reference numbers in each grid square; 
the forest ellipse fitting parameters. 

Each part is created, respectively, by one of the four util- 
ity programs: Map Maker, Hill Maker, Woods Maker, and Forest 
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Maker. When these parts are being created on the Apple 
Disk, the same general file name must be used to identify 
the data's terrain (eg. FOLDA) . Each of the utility pro- 
grams will add its own suffix to the general file name. For 
example, for the Fulda terrain, the four utility programs 
will produce files named: 

FULDA. MAP 
FOLDA. HILLS 
FOLDA. WOODS 
FOLDA. FORESTS 

This suffix identifies the file as a particular portion of 
the terrain data, and will allow the Apple computer to 
recognize it from the other files with the same terrain name 
when reading the file for the combat model. 

1. Ma£ Maker Prog ram 

This program will read and write data files contain- 
ing a catalog by grid square of the hills which influence 
the terrain in that grid square. This program produces an 
equivalent file to that in STAR called LIST. H. Map Maker 
produces a Sequential Text File containing: 

The size of the battlefield (normally 10 by 10) (L,W) 

The lower left grid square in 1000 meter units (LX,LY) 

The base elevation of the battlefield (HL (0,0,0)) 
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For each grid square I,J: 

(1) The number of hills influencing the terrain 
(HL(I,J,0)) 

(2) The reference number of these hills (HL(I,J,K)) 

The user has the option to make a new file or access 
an old file for updating or changes. For new files, the 
program will start with the lower left grid square and pro- 
ceed from west to east then south to north. The data can be 
saved after the information for a grid square has been 
entered. For old files, the program will access it, print 
out the file on the screen and ask for the starting grid 
square. The grid square index is used to identify in which 
grid square to start. The grid square index is an integer 
from 1 to L in the east-west direction, and 1 to W in the 
south-north direction of an L by W sized battlefield. The 
lower left grid square is 1,1; the lower right is L, 1 ; the 
upper left is 1,W; and the upper right is L,W. 

a copy of the program is attached as Listing D in 
Appendix F and a sample file is attached as Figure 15. 
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2. H ill Maker Pro gr am 



This program reads/writes data files containing the 
fitting parameters for each hill on the battlefield. Hill 
Maker produces a Random Access Text File containing: 

The number of hills on the battlefield (N) 

For each hill I: 

(1) X coord of the center location of the hill (H(I,1)) 

(2) T coord of the center location of the hill (H(I,1)) 

(3) The elevation of the hill top measured from 
zero=sea level (H(I,3)) 

(4) The orientation angle of the ellipse measured in 
degrees counter-clockwise from east to the major axis. 
(H(I,4)) 

(5) The eccentricity, defined as the ratio of major 
axis length to minor axis length (H(I,5)) 

(6) The spread, defined as the distance in meters meas- 
ured along the major axis from hill center to the contour 
line which is 50 meters down from the peak (H(I,6)) 

(7) The maximum height of the normal curve describing 
this hill mass (H(I,7)) 

(8) The vertical distance measured from the peak beyond 
which this hill is not considered in computations (H(I,8)) 
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A copy of the program is included as Listing E in 
Appendix ?, and a sample file is attached as Figure 16. 

3. Woods Ma k er Pr ogram 

This program will read and write data files contain- 
ing a catalog by grid square of the woods which influence 
the terrain in that grid square. This program produces an 
equivalent file to that of Map Maker. Woods Maker produces 
a Sequential Text File containing: 

The size of the battlefield (normally 10 by 10) (L,W) 

The lower left grid square in 1000 meter units (LX^LY) 

For each grid square I,J: 

(1) The number of forests influencing the terrain 
(PL(I,J,0)) 

(2) The reference number of these forests (FL(I,J,K)) 

Woods Maker contains the same options as those found 

in Map Maker. A copy of the program is attached as Listing 
F in Appendix F, and a sample file is attached as Figure 17. 

Forest Maker P rogra m 

This program reads/writes data files containing the 
fitting parameters for each forest on the battlefield. For- 
est Maker produces a Random Access Text File containing: 

The number of forests on the battlefield (N) 

For each forest I: 
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(1) X coord of the center location of the forest 

(2) Y coord of the center location of the forest 
(F (Irl)) 

(3) The height of the trees in the forest (F(I,3)) 

(4) The orientation angle of the ellipse describing rhe 
forest measured in degrees counter-clockwise from east to 
the major axis (F(I,4)) 

(5) The length of the semi-major axis of the ellipse 
(F(I,5)) 

(6) The length of the semi-minor axis of the ellipse 
(F(I,6)) 

Forest Maker contains the same options as those 
found in Hill Maker. A copy of the program is included as 
Listing G in Appendix F, and a sample file is attached as 
Figure 18. 

B. THE TERRAIN MODEL PROGRAM 

A copy of the model which explicitly plays terrain is 
included in Appendix F as Listing H. A line drawn along the 
edge of the program delineates the terrain model portion 
which was added to the basic model. The additional varia- 
bles required to run the terrain model are also included at 
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Figure 17: Sample File from Woods daJcer Program 
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the end of this appendix. Professor James K, Hartman [Ref. 
14] provides the details of the methodology for line of 
sight and elevation using the STAR terrain. 

equivalent STAR variables are in ( ) 

A1 = computing parameter for a hill (A) 

B1 = computing parameter for a hill (B) 

CR = critical value parameter for a hill (CRIT.H) 

Cl = conversion of orientation angle to radians 
C2 = (spread of a hill) squared 

FI = value of the elevation of a point due to a hill (FI) 
H(I,J) = the J-th parameter for hill I; the parameters are 

1 = x-coord of the center of the hill (XC.H) 

2 = y-coord of the center of the hill (YC.H) 

3 = peak elevation of the hill (PEAK.H) 

4 = orientation angle of the hill (ANG. H) 

5 = eccentricity of the hill (ECC.H) 

6 = spread of the hill (SPRD.H) 

7 = max height of the hill (HT.H) 

8 = cut height of the hill (CDT.H) 

H(I,0) = status of hill I; 

0 = hill I has not been checked for LOS 

1 = hill I has been checked for LOS 
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HL(I,J,K) = the number of the K-th hill in grid square I,J 
which influences the terrain in that grid 
square I = 1 to L; J = 1 to H; K = 1 
to HL(I,J,0) 

HL(I,J,0) = the total number of hills in grid square I,J 

which influence the terrain in that grid square 
HL (0,0,0) = the base elevation of the terrain map 
H2 = temporary variable identifying the hill being checked 
for line of sight 

L = number of grid squares on the terrain map in the 
K direction (NGRIDX) 

LX = x-coord of lower left grid square of terrain map in 
1000 meter units (X.LO.BDEY) 

LX = y-coord of lower left grid square of terrain map in 
1000 meter units (Y.LO.BDRY) 

M$ = name of files containing terrain data 

NH = total number of hills on the terrain map (NHILLS) 

PI = computing parameter for a hill (PXX.H) 

P2 = computing parameter for a hill (PXY.H) 

P3 = computing parameter for a hill (PYY.H) 

QI = quadratic function of the ellipse represnring the 

hill (QI) 
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H = number of grid squares on tbe terrain map in the X 
direction (NGRIDY) 

X = current x-coord of unit whose elevation is being 

computed (in 10(X meter units) (X) 

XB(I) = current x-coord for BLUE unit I 
XS(I) = current x-coord for RED unit I 
XS = distance in x direction between unit location and 
center of hill mass (XS) 

XI = relative number of grid squares in x direction of a 

BLUE unit location with respect to the lower left 
grid square 

X2 = relative number of grid squares in x direction of a 
RED unit location with respect to the lower left 
grid square 

X3 = number of grid squares in x direction between a BLUE 
unit and a RED unit 

X4,T4 = current grid square being checked for line of 
sight in LOS subroutine 

Y = current y coord of unit whose elevation is being 
computed (in 100 meter units) (Y) 

YB (I) = current y-coord for BLUE unit I 
YR (I) = current y-coord for RED unit I 
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YS = distance in y direction between unit location and 
center of hill mass (YS) 

Y1 = relative number of grid squares in y direction of a 
BLOB unit location with respect to the lower left 
grid square 

Y2 = relative number of grid squares in y direction of a 
RED unit location with respect to the lower left 
grid square 

Y3 = number of grid squares in y direction between a BLUE 
unit and a RED unit 

Z = temporary variable for the elevation at a point on a 
hill (Z) 

ZB (I) = elevation of BLOE unit I 

ZR (J) = elevation of RED unit J 
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APPENDIX E 

SYNOPSIS OF THE POSTPR O CESSOR 
The following appendix is a sunnary of the work done by 
other students as a part of a joint project for OA 4655. 

The authors gratefully acknowledge the work done by Captains 
Anbrose R. Hock and Steven L. Maddox in developing this 
postprocessor to reduce the anount of work required in 
extracting the necessary data fron STAR. Their work is 
included here in this appendix since it is an integral part 
of the parameter estimation process. 

The postprocessor developed for this thesis is a Sim- 
script language program that provides to the user a summary 
of the data generated from the Simulation of Tactical Alter- 
native Responses (STAR) Model. Additionally, the postproces- 
sor performs some routines that assist the user in the data 
analysis of the model. The STAR Model is a Simscript lan- 
guage program that simulates combat between two combined 
arms teams in a combat environment that includes field 
artillery fire and electronic warfare. 

Since the model assumed in the project only required 
data generated by a tank vs tank battle, it was necessary to 
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change certain aspects of the STAB model. Thus it was 
necessary to "turn off" the field artillery and Electronic 
warfare modules within the STAR model. Also the non-tank 
units needed to be deleted. Referring to figure 19, the 
Electonic Warfare code(EWCODE) was deleted through the use 
of a comment card, and a Electonic Warfare off(EWOFF) 
subroutine was inserted. 
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Figure 19: Electronic Warfare Deletion 

Thus, whenever the main program called for the usage of EW, 
the EWOFF subroutine automatically returned the program back 
to the point where the EW routine was called and the simula- 
tion continued with no Electronic Warfare generated. The 
field artillery within the STAR model was cancelled by 
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having the simulation create FA fire at time equal to 5000 
time units (see Figure 20). 

Schedule a FABEGIN in 5000.0001 time units 

Figure 20: Field Artillery Deletion 



Since the simulation was scheduled to run only 2500 time 
units, no field artillery was generated. Finally, the 
internal logic of the STAR model was used to "destroy” the 
non-tank units within the simulation by changing their basic 
load to zero (see Figure 21). This action created situa- 
tions, where the tank elements did not consider the non-tank 
elements as lethal or dangerous targets and thus were 
ignored. 



1 


1 


1 


1 


0 


38 


1 


1 


0 


17 


1 


1 


0 


100 


0 


0 


0 


0 


2 


3 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


2 


4 


1 


0 


1 


10 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


3 


6 


0 


0 


0 


0 


1 


0 


1 


2 


0 


0 


0 


0 


0 


0 


0 


0 


4 


1 


50 


20 


0 


50 


60 


30 


0 


50 


0 


0 


0 


0 


0 


0 


0 


0 


1 


7 


1 


1 


0 


22 


1 


1 


0 


18 


1 


1 


0 


100 


0 


0 


0 


0 


2 


8 


1 


0 


1 


4 


0 


0 


0 


0 


1 


1 


0 


1 1 


0 


0 


0 


0 


1 


1 


1 


1 


0 


38 


1 


1 


0 


17 


1 


1 


0 


100 


0 


0 


0 


0 


2 


3 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


2 


4 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


3 


6 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


4 


1 


50 


20 


0 


50 


60 


30 


0 


50 


0 


0 


0 


0 


0 


0 


0 


0 


1 


7 


1 


1 


0 


22 


1 


1 


0 


18 


1 


1 


0 


100 


0 


0 


0 


0 


2 


8 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 



Figure 21: Alteration of Basic Load 



Output from the STAR model can be divided into three 
categories: 
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Opiate list of attacJcer/defender status. 

Shot data 

Significant event listing (i.e, movemenr, EH and ?k 
events) 

Only the shot data from the STAR model is needed to provide 
the necessay data for computations. Figure 22 is a sample 
of this output. In order to utilize this data, the internal 
output control of the model was used to place the shot data 
in a mass storage location within the computer to be later 
manipulated by the postprocessor. 

The flow chart for the postprocessor is contained in 
Figure 23; a copy of the program itself can be found in 
Appendix F as Listing I. The program flow is very basic and 
needs little explanation. From the mass storage location of 
the computer, the postprocessor reads pertinent information 
into the SHOTLIST. Once all of the data has been read, the 
program processes the data into a casualty list (CASLIST) 
keying from a status of ’’dead” from the Shotlist. Figure 24 
is a sample of the casualty list output. Once CASLIST is 
created the postprocessor performs histogram computations 
and graphs. The postprocessor can print two types of 
histograms. Figure 25 contains the Simcsript Library pro- 
gram for the histogram. Additionally, the HISTG routine in 
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Figure 22; STAR Output Shotlist 
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FORTRAN can be accessed to produce a histogram of the data. 
After the histogram data is generated, the postprocessor 
calculates the attrition rate coefficients using the formu- 
las discussed in Section III-E of this thesis. 

Bri e f descript ion of the postprocessor (see this appen- 
dix for variable description) : 



Preamble 


• 

• 




lines 


1 - 9 


definition of permanent entities. 


lines 


10 - 39 


definition of temporary entities. 


lines 


40 - 54 


definition of data to be collected for 


histograms. 




Main ; 






line 


2 


variable definition. 


line 


3-4 


system parameters read. 


lines 


13 - 31 


creation of shotlist, reading data from 


file. 






lines 


43 - 89 


creation of casualty list, file '’dead”. 


lines 


90 - 125 


processing of histogram data 


lines 


125 - 186 


processing of attrition rate coeffi- 


cients 


• 
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SHOT LIST DATA FOR CASUALTIES FOR THIS STAR SIMULATION 



NUM 


TIME 


TARGET SYS 


-WPN 


TYPE TIME BETWEEN CASUALTIES 


1 


322 


59 


1-7 


0 


2 


334 


46 


1-7 


12 


3 


46 8 


58 


1-7 


134 


4 


57 4 


28 


1-1 


106 


5 


581 


42 


1-7 


7 


6 


591 


27 


1-1 


10 


107 


1210 


72 


1-7 


6 


108 


1238 


69 


1-7 


28 


109 


126 5 


73 


1-7 


27 


1 10 


130 3 


75 


1-7 


38 


111 


1323 


19 


1-1 


20 


112 


1421 


76 


1-7 


98 






Figure 


24: 


Casualty List 



131 



f 



NOMBER 
BLMO = 
N0M3ER 
RDMU = 



AHRT = 



BHAT = 



BLHST 



3LMU = 



OF BLOE KILLED = 28 

47,25000 BLSIG = 299.259 
OF RED KILLED = 84 

13,24096 RDSIG = 219.106 



HISTOGRAM OF BLOE AND RED CASOLTIES 
(DATA BASED ON TIME BETWEEN BLOE CASOLTIES 
AND TIME BETWEEN RED CASOLTIES) 



INTERVAL=10 ONITS 


#BL0E CAS 


#RED CAS 


1 


3 


11 


2 


4 


8 


3 


4 


20 


4 


0 


4 


5 


0 


8 


6 


0 


4 


7 


0 


0 


8 


0 


8 


9 


0 


0 


10 


0 


0 


11 


0 


4 


12 


4 


8 


13 


0 


4 


14 


4 


4 


15 


0 


0 


16 


0 


0 


17 


0 


0 


18 


0 


0 


19 


4 


0 


20 


0 


0 


21 


5 


0 



Figure 25: simscript Histogram 



Li^ of Variable s For Postprocessor 
estimate of the attrition coefficient for RED 
attriting BLOE 

estimate of the attrition coefficient for BLOE 
attriting RED 

Simscript routine for the histogram of BLOE time 
between casualties 

Simscript routine for calculating the mean time 



132 



1 



between casualties 



BLSHTHST = Simscript routine for the histogram of BLUE 
shot times 

BLSHTIG = Simscript routine for calculation of the 
standard deviation of BLUE shot times 
BLSHTMU = Simscript routine for calculation of the mean 
of BLUE shot times 

BLSIG = Simscript routine for calculation of the 

standard deviation of BLUE time between casualties 
BLUE. SHOTS = time of the previous BLUE shot 
BLUE. TIME = time of the previous BLUE casualty 
CAS = casualty 

CASHST = Simscript routine for the histogram of 
toral time between casualties 
CASIG = Simscript routine for calculation of the 

standard deviation of the total time between 
casualties 

CASLIST = array containing list of casualties 
CASMU = Simscript routine for calcualtion of the mean 
for the total time between casualties 
CBLUE = number of surviving BLUE force elements 
CK = integer used in the calculation of the attrition 
coeff icient 
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1 = BLUE casualty 
0 = RED casualty 

CRED = number of surviving RED force elements 

CURTIME = time of previous casualty 

FNAME = firer*s name 

FSWTYPE = firer’s system-weapon type 

J.BL = number of BLUE casualties (used in BLHST routine) 
J.RD = number of RED casualties (used in RDHST routine) 
JBLHST = number of BLUE shots(used in BLSHTHST routine) 
JCAS = number of casualties (used in CASHST routine) 
JRDSHT = number of RED shots (used in RDSHTHST routine) 
JSHOT = number of shots( used in SHOTHST routine) 

JTOT = number of casualties (used in TOTHST routine) 

LIM = array of time periods that represent phases of 
the battle 

LIMA = lower time limit of phase 
LIMB = upper limit of phase 
LP = partial derivative of p 
LQ = partial derivative of q 
MAXBLUE = initial BLUE force size 
MAXRED = initial RED force size 
MINP = minimum of p 
MINQ = minimum of q 
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MK1 = number of RED elements in time interval 

NK1 = number of BLUE elements in the .time interval 

NLIM = number of time intervals to be evaluated 

NSHOT = number of shots in the total battle 

NTABLE = number of tables (battles) to be evaluated 

NUM = counter fo number of casualties 

NUHBL.CAS = number of BLUE casualties 

NUMRD.CAS = number of RED casualties 

NUMCAS = casualty entity number 

P = probability of not acgiring a target 

PHAT = estimated p 

PROVAL = process variable 

0 = stop 

1 = process 
QHAT = estimate of q 

RANGE = range of firer to target 

RDHST = Simscript routine for the histogram of RED 
casualt ies 

RDMU = Simscript routine for calulation of the average 
time between RED casualties 

RDSHTHST = Simscript routine for the histgram of total 
red shots 

RDSHTMU = Simscripr routine for calculation of the average 
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red shot time 



SDSHTSIG = Simscript routine for the calculation of the 
standard deviation of the RED shot time 
RDSIG = Simscript routine for calculation of the 

standard deviation for RED time between casualties 
RED. SHOTS = time of previous RED shot 
RED. TIME = time of previous RED casualty 
RONNO = run number 
SHOT = shot 

SHOTHST = Simscript routine for total shots fired 
SHOTLIST = array of shot times 

SHOTMU = Simscript routine for calculation of the average 
time of shots 

SHOTSIG = Simscript routine for calculation of the 
standard deviation of total shots 
STATUS = status of shot 
STATUS. CAS = status of casualty 
T.BTWN.CAS = time between casualty 
T.CAS = casualty number 

TBLUE = number of BLUE casualties within the interval 
TGT.CAS = target casualty 
TIME = time of shot 
TIMECAS = time of casualty 
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TIMELIM = maxtime of battle when a casualty occurred 



TMCAS = time of previous casualty 
TMRDCAS = time of previous RED casualty 
TOTAL. CASTIME = total casualty time 
TOTAL. SHOT. TIME = total shot time 

TOTAL. TIME = total Casualty time for the creation of 





T.BTWN.CAS 


TOTHST = 


= Simscript routine for the histrogram of 



total casualty time 



TOT MO = 


Simscript routine for calculation of the 



average casualty time 



TOTSIG = 


: Simscript routine for calculation of the 



standard deviation of the total casualty time 



TRED = 


:otal RED casualties within the time interval 



TSOB = time counter 



TSWTYPE 


= target system-weapon type 


TTBLOE = 


= total BLUE casualties for XX battles 


ITRED = 


total RED casualties for XX battles 


TTCAS = 


total casualties for XX battles 



TTT = number of shot 

TYPE. CAS = system weapon type of the casualty 
VAR = attrition coefficient variables 
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APP SNDIX F 

LISTI NG OF C^PajSS PH OGPAaS 
A. DAMSTAC MODEL PROGRAM 



i0 


REH 


UPDATED 24 FEB 82 






20 


REH 


DETERMINISTIC AC-ieRE6ATE0 MODEL OF STQR 


HN 


THE 


21 


REM 


APPLE COMPUTER < DAMSTAC) 






30 


DIM 


XB< 5 ) /r'B< 5 > .KR< 5 >, VR< 5 5^10) »3V< 5 , 


1 0 > 


,RX<5,10 


40 


DIM 


BPC 5 ) .RP< 5 > ,SB< 5 ) ,SR^ 5 ) ,RVr 5 . 10 






50 


DIM 


SF< 5 > ,RF< 5 > ,8T< 5 ,5 ) ,RT< 5,5 > ,nB< 5 -• ,0R< 5 > 






60 


DIM 


BR<5>,RR<5) 






70 


DIM 


8lC5>,ftl< 5> 






80 


T s 0;OT s 30; MT = 1500 







80 Ds = 


CHRS < 4 > 


100 


P s 0 


.38:Q s 0.33 


110 


PI - 


3. 14 1582654: PR = 0 


120 


A = 2 


:.4;B = 1.2 


130 


BE * 


0,3: RE * 0.3 


140 


BM = 


,5;RM = ,5 


150 


KK = 


0 


160 


Nk‘ = 


MT / DT 


ISI 


REM 


22=FLA6 FOR TVPE OUTPUT 


162 


REM 


22-1 0IUE3 UNIT 


183 


REM 


LOCATIONS •$< SI2ES 


164 


REM 


22=0 6IUES TIME TOTAL 


165 


REM 


RED :?< SLUE FORCE LEUELS 


166 




0 


170 


OIM 


B0< NK > .R0< NK > 


130 


HOME 


181 


REM 


READ DATA FROM FILE FOR 


182 


REM 


BLUE FORCE. 


133 


REM 


DATA INCLUDES: 


134 


REM 


MssNUMBER OF UNITS 


135 


REM 


01=NUMBER OF COORD PTS 


186 


REM 


BX,8V=6RID OF COORD PT 


13? 


REM 





180 IHP«JT “NPHE OF FILE '=^OR RLME OOTp-7- '• :NS 
800 PRINT OS; "OPEN ";NS;".eUJE” 

8 10 PR I NT OS ; " READ '* ; NS ; •• . BLUE '* 

220 INPUT M 

230 FTiR I = 1 TO M 

240 INPUT 81 

250 FOR ..! = I TO 81 

260 INPUT B:2< ! .J>; INPUT BV-: I , J 

270 NEI'^IT : ^^6’XT 

280 PRINT OS; “CLOSE NS;" , 9LUE’‘ 
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ZSl ftSH SET IN ITT PL UHTT 
282 REM LOCPTTON EQUPL TO 
233 REM 1ST COORO PT. 

284 REM 

238 FOR T = t TO M::2e<I> = 8r<( I ^ T >: VBC I > s 8VC I , t ■> s PP< I > = 1 
308 PRINT 81: FOR JJ = 1 TO Bl: PRINT 8X<' I ^JJ '• »ev-r I >: NEXT JJ 

310 PRINT “SIZE OF BLUE UNIT “?T;; INPUT SB< I '- 

320 81< I > = S8< I j 

321 REM 

322 REM SUM INITIPL FORCE SIZE 

323 REM 

330 80< 0 > = 60< 0 > *■ SB< T > 

331 REM 

332 REM RPTE OF BLUE MOUEMENT 

333 REM EQUPLS 8 

334 REM 

340 8R< I > = 0 

350 NEXT 

360 PRINT 

361 REM 

362 REM REPCi OPTP FROM FILE FOR 

363 REM RED FORCE. 

364 REM SPME FORMAT PS BLUE 

365 REM 

370 INPUT “NPME OF FILE FOR RED OPTP'? “ ;N« 

380 PRINT 0:5; “OPEN “;N«;“.RED“ 

390 PRINT OS; “REPO “;NS;“.RED“ 

400 INPUT N 

410 FOR I = 1 TO N 



420 INPUT R1 

J-30 FOR .1 = 1 TO R1 

440 INPUT RXa.J>; INPUT RV< I,J> 

450 N0<T ; NEXT 

460 PRINT D:S; “CLOSE “;N5:'‘.RED" 

470 FOR I = I TO N:XR<I> = RX< I , r): VP< I > - R**'< I . t RP< r> = 1 
480 PRINT “SIZE OF RED UNIT INPUT SR< I> 

•*90 Rl( I > = SP/ I 'f 

500 R0(’ 0 > =s R0< 0 > + SR'' I > 

510 RR< I > = 5 
520 NEXT I 

530 PRINT ; input “DO VOU NPNT THE PRIN'rnp v/n :• “ :P5 

535 IF Pi = “V" THEN PR = 1. 

540 rHjME ; UTPB 3 

550 IF PR = 1, THEN PR# 1 

560 PRINT 

564 REM 

565 REM PRINT INITIPL OPTO 

566 REM 

570 PRINT “TIME";; POKE 36.7; PRINT “PU..1F”;: POKE 36.20; PRINT “RED" 
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580 PR# 0 

584 REM 

585 REM WHICH OUTP«JT TVPE? 

586 REM 

580 IF 22 = 1 THEN t?60! 8CiTQ P10 

600 iWSUe 1490 

610 REM MHIN PROeRHH 

620 T = T + OT 

624 REM 

625 REM CHECK IF MAX TIME 

626 REM IS EXCEEDED 

627 REM 

630 IF T < = MT TH£h< 840 

635 PRINT -END OF 80TTLE DME TO TIME"; 603UB 1540: 
640 60SUB 790: REM MOUEHENT SUBROUTINE 
650 I30SUB 1170: REM ATTRITION SUBROUTINE 

654 REM 

655 REM INCREMENT THE TIME 
856 REH PERIOD INDEX 

657 REM 

660 K.K = KK + 1 

670 60<KK> = 0:R0<KK :' = 0 

674 REM 

675 REM SUM BLUE FORCE LEUEL 

676 REM FOR THIS TIME PERIOD 

677 REM 

530 FOR I = I TO M 



690 


60< KK 


= 80< KK > + SB<' I 


?00 


NEXT 


I 


704 


REH 




705 


REH 


SUM RED FORCE LEUEL 


706 


REM 


FOR THIS TIME PERIOD 


707 


REM 




710 


FOR 


J = 1 TO N 


720 


R0< kk; 


> = R0C KK > SRc J 


730 


NEXT 


J 


734 


REH 




735 


REM 


WHICH TYPE OF OUTPUT'- 


736 


REH 




740 


IF 2 


2 = 1 THEN GOSJJB 1760 


750 


GOSUB 1400 



END 
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752 REH 

753 R£M CHECK IF EHO OF 8HTT1 F 

754 REM - CR I TER IP FOR FORCE 

755 REM LEOELS IS EXCEEDED 

756 REM IF CR I TER IP IS EXCEEDED 

757 REM 30T0 SPUE RESULTS 

758 REM 

768 IF e8<KK:) > BE 4 88<8> THEM 778 

765 PRINT : PRINT "EHO OF BPTTLE DUE TO PTTRITION OF BLUE" 

766 0OSUB 1540; END 

770 IF R0<KK) > RE 4 R0<0> THEN 788 

775 PRINT : PRINT "END OF BPTTLE DUE TO PTTRITION OF RED" 

776 60SUB 1540; END 
730 iKiTO 620 

790 REM MOUEMEh^T SUBROUTINE 

794 REM 

795 REM FOR RED UNITS 

796 REM 

300 FOR I s 1 TO N 

304 REM 

305 REM SET LPST COORD PT 



306 


REM 




310 


J * RP<I > 


313 


REM 




314 


REM 


COMPUTE COMPONENTS OF 


315 


REM 


OISTPNCE BETWEEN NEXT 


316 


REM 


COORD PT PND CURRENT 


317 


REM 


LOCPTION 


313 


REM 




320 


DV s 


RV< I , J 4 1 > - VR< I > 



330 OX * RX<; I , J + 1 ) - XR< I > 

333 REM 

334 REM DO SEOMETRV TO SET 

335 REM CORRECT PN6LE (PN^ 

336 REM 

340 IF OX ~ 0 PNO OV >• 0 then PN = PT ••• SOTO 330 
350 IF ox = 0 PNO OV < 0 THEN PN = 3 4 PT ' i50TO 330 
360 PN = PIN <DV / DX ) 

370 IF OX < 0 THEN pN = PI + PN 
373 REM 

^74 REM COMPUTE OH, OISTPNCE 

375 REM TO BE TRPUELLED IN 

376 REM THIS TIME PERIOD 

377 REM 

330 01 = OT 4 RR< I ■) * COS PN 

390 02 = OT 4 RR< I > 4 SIN (PN;:* 

300 OH = SQR <01 2 4 02 2> 



141 



804 


REM 






805 


REM 


COMPUTE 01 » 01 


STANCE 


906 


REM 


TO NEXT COORD 


PT 


807 


R£M 


FROM CURP^T location 


808 


REM 






810 


01 s 


3QR <0X ^ 2 


OV •' 2> 


811 


REM 






812 


REM 


IF 01 < OH, MC 


:iUE UNIT 


813 


REM 


ONLV THE AMOUNT TO 


814 


REM 


NEXT COORD PT 




815 


REM 


IF 01 IS USED. 


, UPDATE 


816 


REM 


LAST COORD PT 




817 


REM 






820 


IF 01 > =* OH THEN 


950 


830 


01 s 


OX; 02 = OV 




840 


RP< I 


> = J -f 1 




850 


VR< I 


■> = TNT (VRCT :- 


02 ':' 


860 


XR< I 


> = INT(XR<T> 


+• 01 '> 


870 


NEXT 




874 


REM 






875 


REM 


00 THE SAME FOR BLUE 


876 


REM 






880 


FOR 


I * 1 TO M 




885 


IF 


BR( I > = 0 THEN 


1150 


890 


J = 


8P< I > 





Ifcjgie DV = 8V< I,J + I > - ve< I :■ 

1010 OX - 8)<c I,J + r> - :<8< T> 

1020 IF OX = 0 mNO OV > 0 THEN ON = PT / i: r,riTn 1,0150 

1030 IF OX s 0 AND OV <’ 0 THEN AN = 3 * PT / 2 : 0OTO 1060 

1040 AN = ATN <DV / OX > 

1050* IF OX < 0 THEN AN = PT + AN 
1060 01 = OT # 8R< I > * COS <AN > 

1070 02 = OT * eRCI -' ^ SIN cAN^ 

1080 OH = 3QR 'roi " 2 + 02 2) 

1090 01 = 3QR ' OX •' 2 + OV 2> 

1100 IF 01 > = OH THEN 1130 

1110 01 = OX; 02 = OV 
1 120 9P< I > - .J -t- 1 
1130 VBt:i> = I NT 'rVB'ri) + 02:- 
1 i-AO xe< I = I NT < XB< I 01 :> 

1 150 NEXT 
1160 RETURN 
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1170 


REH 


RTTRITTON SUeROUTINE 


1172 


REH 




1173 


REM 


ZERO OUT QTTRITTON 


1174 


REH 


UHLUES FOR RLL UNITS 


1175 


REH 




1160 


FOR 


I = 1 TO MieF<I> = 0:Oe< 


1 1S0 


FOR 


J = 1 TO N:RF<J) * 0:DR<. 


1132 


REH 




1133 


REM 


CHECK RRN6E BETWEEN 


1194 


REH 


OPPONENTS, IF RANGE 


1135 


REM 


<=: 2700 ATTRITION WILL 


1196 


REM 


OCCUR. REDUCE RATE OF 


1137 


REH 


RED MOUEMENT TO 2 H-^' 


1138 


REH 




1200 


FOR 


I = 1 TO M 


1210 


FOR 


J = 1 TO N 


1220 


8T( I 


,J> = 0:RT<J,n s 0 


1230 


R6 = 


3QR < < Xe< I > - XR< J :> ■> ^ 



= NEXT 
s ^3; NEXT 



+ •' I - 



VP< J ■> > 



1240 IF Ri3 > 2700 THEN 1310 
1250 RR(.J> s 2 



1251 REN 

1252 REH 

1253 REM 

1254 REM 

1255 REM 

1256 REM 

1257 REM 



INCREMENT THE NUMeER 
OF TOTS FIRED RT 8V 
EACH UNIT. COMPUTE 
HTTP I TI ON CfiUSED 8V 
RED ON BLUE S< UICE 
UERSR 



1258 


REH 












1270 


8F< I > = BF< I > 




1 








1280 


RF<J) = RF<J> 


+> 


1 








1235 


IF RR<.J) = - 6 


1 THEN BKI 


,J :• 


— 


0: 'R»*'TO 1295 


1290 


8T<I,J> = A / 


N 


* (■ 1 - P > 




SBC 


I>) 4 SRC-r- 


1235 


IF 9R<r> = 3 


THEN RT^J,!') 




0: 


ROTO 1310 


1300 


RT<J.I> = 8 


H 


4(1 - ( Q > 




SRC 


.J y ') •* SRC I > 


1310 


NEXT : NEXT 












1311 


REM 












1312 


REM SUM TOTAL 


ATTRITION 








1313 


REM FOR EACH I. 


IN IT 








1314 


REM 












1320 


FOR I = 1 TO 


M 










1330 


FOR J a 1 TO 


N 










1340 


IF RF':.J> a 0 


OR 3F< I > a 0 


THPN 


t370 


1350 


08C I ) a DB< I 




8T( I , J •) / 


RF( 


. J ■> 




1 360 


opxj:> a oR< j> 




RT<'J,D ••• 


3F<’ 


I 




1370 


NEXT : NEXT 
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1371 

1372 

1373 

1374 

1375 

1376 

1377 
1373 
1379 
1330 



ftHH 

3EH COHPUTE NEW IINTT 
ftEH STRENGTH. SUBTRACT 
KEM unit STTRITTfiW coriH 
REH CURRENT UNIT STRENGTH 
REH IF UNIT STRENGTH T-fi'iRq 
REH NEGGTIUE> MPKE IT 
REH EQUPL 0. 

REH 

FOR I s 1 TO M 



1390 S8< I > = SBC I > - OBC I > 

1400 IF SBC T > < 0 THEN SR'' T ') = 

1410 IF SBC I > < = BM # BIC I > THEN BRC I 



1 h20 next 



1430 FOR J = I TO N 






1440 3RCJ> = 3RCJ> - ORCJ > 

1450 IF SRCJ ) C 0 THEN SRCJ > = 0 

1460 IF SRCJ ) < = RH 4 PI' I - THFW PP<J > = - 3 

1 470 NEXT 

1460 RETURN 

1490 REH PRTNTOUT--'SUHMGRV 

1500 IF PR = 1 THEN PR4> 1 

1510 PRINT T;: POKE 36.7? PRINT pp^ckk-.t 

1515 POKE 36,20: PRINT R0CKK) 

1520 PR# 0 
1530 RETURN 

1540 REH SPUE RESULTS IN P FTl.E 
1550 PRINT 

1555 INPUT “DO VOU WPNT P FILE '*'F twtc pun^c nqTp'? 
1560 IF PS = “N" THEN 1750 
1570 OS = CHRS <4> 

1530 INPUT “NPME OF FILE FOR FORCE i FiiCL DPTP-*' ” ;FS 
1590 PRINT OSi“OPEN '' jFS; “ . COEFS’* 

1600 PRINT 0Si“WRITE “ ;FS;‘' . COEFS" 

1610 PRINT P: PRINT 8; PRINT P: PRINT Q 
1620 PRINT OS;“CLOSE ” :FS; “ . COEFS" 

1630 PRINT OS; “OPEN “;FS:“, RESULTS" 

1640 PRINT OS;“WRITE “ ;FS;“, RESULTS" 

1650 PRINT NK + 1; PRINT NK + 1 

1660 FOR I = 0 TO NK 

1670 T = I * .3 

1630 PRINT T: PRINT 60C I > 

1690 NEXT I 

17061 FOR 1=0 TO NK 

1710 T = T ^ 3 

1720 PRINT T: PRINT R0< I'-* 

1730 N0<:T I 

1740 PRINT OS;"i::LOSE “ ;FS; RESULTS" 

1750 RETURN 



“ :PS 
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1 



i76w 

1770 

1730 

1790 

1800 

1310 

1320 

1330 

1840 

1350 

1830 

1870 

1330 

1330 

1300 

1310 

1320 

1330 

1340 



REH PR I NTOUTx'SUHHPR V 
HOME : UTAB 3 
IF PR s 1 THEN PR# I 
PRINT "TIME IS “;T 
PRINT : PRINT 
PRINT "BLUE FORCES" 

PRINT "UNIT"; TOBC 7 >;"><:cn0R0" ; TCiB< t3>;“WnRn“; Tqe< 23 >; “SIZE" 
FOR I = 1 TO M 

PRINT " “;I; TPB< 7>;MB<I>; TQB*' 15>;Ve<I:’: tqb < 24:.;S8d> 

NEXT 

PRINT : PRINT "RED FORCES" 

PRINT "UNIT"; Tfie< 7 > ; "XCOORD" ; TflR^ t=i i-VCOORn" ; TCiB< 23 >; "SIZE" 
FOR I = 1 TO N 

PRINT “ ";I; Ti38< 7>;XR<I>; TQB< t5':';VR<r>; TQ8< 24;;-;SR<I> 

NEXT 

PRINT ; PRINT 
PR# 0 

INPUT "HIT RETURN TO CONTINUE” ;Zf 
RETURN 
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B. FOaCB BAKER PROGHAB 



lti REM FORCE OPTQ-FII.E MPkER 
20 DIM BX< 5, 10 ) ,8V< 5, 10 ) ,ei< 10 •• 

•50 GOTO 1 70 
40 HOME 

50 INPUT “NPHE OF FILS'? “ ;NS 

G0 OS = CHRS ' ^ j 

70 PRINT OS; “OPEN “ ;NS;“ . “ ;MS 

30 PRINT OS;“WRITE “;NS;“.“;MS 

30 PRINT M 

100 FOR I = 1 TO M 

110 PRINT ei<i> 

120 FOR J = 1 TO Bl<l> 

130 PRINT 8X‘<I,J>: PRINTeV<t,J^ 

140 NEXT .J 

150 NEXT I 

160 PRINT OSj“CLOSE “;NS;“.“;MS 

170 HOME 

180 UTPe 5 

190 PRINT “MENU:" 

200 PRINT “1 MPKE FILE FOR BLUE FORCES" 

210 PRINT "2 MPKE FILE FOR RED FORCES" 

220 PRINT "0 QUIT" 

230 PRINT ; INPUT "WHICH? ";Z 

240 IF 2 s 0 THEN END 

250 IF 2 = I THEN MS = “BLUE" 

2S0 IF 2 = 2 THEN MS = “RED" 

270 HOME 

280 PRINT "WHPT IS THE FORCE -SIZE OF ";MS;i INPUT “ "jM 

290 FOR I = I TO M 

•300 PRINT ; PRINT "FOR UNIT ":I 

310 PRINT “HOW MPNV COORD I NOTING POINTS “ ; 

•515 INPUT "(INCLUDING STORT LOCHTTONS :•? " ;et( r> 

320 FOR ...I = I TO BK; !':• 

330 PRINT 

340 PRINT "ENTER GRID COORDS FOR COORD T NOTING P-OINT * “ 
350 INPUT “X-COORD IS " jBXC I -■ 

360 INPUT "V-^::OORD IS “ ;BV<' T ;■ 

370 NEXT J 

380 NEXT I 
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390 HOHE 

400 PRINT ; PRINT Tfi8< DOTP CHErK 

410 PRINT : PRINT : PRINT 

420 PRINT “UNIT PT“; TPe-^ 9>;“X-H:Cn*iR0“ ; TPB< t?>;“V-H:OORD“ 

430 FOR I = 1 TO M 
440 FOR J = I TO 91< I 

450 PRINT TPe< 2>;I; T|:)B< TOe< TB8< 13>;8V<I,J> 

460 NEXT J 
4?0 NEXT I 

*80 INPUT “PNV CiHPNeeS? <V/N> “;2‘S 
490 IF ZS = “N“ THEN 40 
500 PRINT 

5i0 INPUT “WHICH UNIT^' “;I 
520 INPUT “WHICH POINT'' “;J 
530 INPUT “X-«:4D0R0 IS “;8Xa,J> 

540 input “V-COORO is “;8V<I,J> 

550 SOTO 390 
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c 



HESOLTS SEADSa PROGRAM 



iti HCiHE ; UTA6 3 
20 — CHRf (. 4 > 

30 INPUT “HaHE OF FILE FOR FORCE '.EUEU DRTP'? ” :Ff 
40 PRINT OS i “OPEN “;FS;“.COEFS“ 

50 PRINT OS; "REPO ";FS;“.COEFS" 

80 INPUT P: INPUT 8: INPUT P: INPUT Q 
T0 PRINT OS; "CLOSE “ ;FS;" , COEFS" 

80 PRINT DS;“OPEN " ;FS;“. RESULTS'* 

30 PRINT DS;"REPD ";FS;“. RESULTS" 

100 INPUT N; INPUT N 

110 NK = N - 1 

120 DIM B0<NK>,R0<NK> 

1 30 FOR I = 0 TO NK 
140 INPUT T: INPUT 80(1 
150 NEXT I 
180 FOR 1=0 TO NK 
170 INPUTT: INPUT K0( I> 

130 NEXT I 

130 PRINT Os;"CLOSE " ;FS;“ . RESULTS" 

200 INPUT "PRINTER ON? “;PS 

210 HOHE : UTPe 3 

220 IF PS = -t'" then PR# 1 

230 PRINT “P = " ;P 

240 PRINT “8 = ";e 

250 PRINT "P =s “ ;P 

260 PRINT "Q = “;Q 

270 PRINT 

230 PRINT "TIME":: POKE 38*7: PRINT “SLUE"; 

235 POKE 38,20: PRINT "RED" 

290 FOR KK = 0 TO NK 
300 T = 30 * kK 

310 PRINT T;: POKE 38,7: PRINT 80<KK:::>; 

315 POKE 36,20: PRINT PfXKK) 

320 NEXT KK 
330 PR# 0 

340 ENO 
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map KAK2B PROGRAM 



10 REM HPiP CflTRLOG MPKER/REfiOEH 
20 DIM HL<10»10»12> 

2.0 HOME ; UTPB 3 
40 PRINT "HENU:“ 

50 PRINT “1 MPP MPKER" 

60 PRINT "2 MhP REPOER" 

73 PRINT "0 QUIT" 

30 PRINT : INPUT “WHICH^ ‘‘;K 

30 IF X = 0 THEN END 

100 ON X SOTO 120 #.300 

110 REM MfiP CPTPLOG MPKER 

120 HOME : UTPB 2 

130 PRINT TP8( li:>;"^HHf MPP MPKER 

140 PRINT ; PRINT 

150 INPUT "NEW OR OLD MPP^ <N/0> " 

160 IF Z$ = "N" THEN II = 1:JJ = 1: SOTO 220 
170 INPUT "NPME OF MPP"? ";M$ 

1 30 DS = CHRS '• 4 ) 

130 SOSUB 620 

200 INPUT "ENTER STPRTTNG SRID 60UPRE 
210 SOTO 250 

.220 INPUT "SIZE OF MRP <1000 M SOUPRES'-'^ ";I.#W 

230 INPUT "ENTER LOWER LEFT SRID SQUPRE <ES. 50#.33> ";LX#L 

240 INPUT "8PSE ELEUPTION OF MPP^ ":HL<0#0#0> 

250 FOR I — II TQ L 
260 FOR J s .JJ TO W 
270 PRINT : PRINT 

230 PRINT "FOR SRID SQUPRE ";I;".";J 

2S0 INPUT "NUMBER OF HILL MPSSES'^ " ;HL< I # J#0 >: KK = HL<I>J# 
•300 IF KK - 0 THEN 3.S0 

310 PRINT "ENTER THE ";KK;“ HILLS FOR hRID POUPRE ";I#"#"; 
320 FOR K = 1 TO KK 

330 PR I NT "HILL " ; K ; " IS " ; ; I NPUT ml.:; i , j, k :> 

340 NEXT ;< 

■350 JJ = 1 

360 PRINT ; INPUT "SPUEP CV/N'VjPS 
370 IF = "V" THEN 400 
■380 NEXT J 
330 NEXT r. 



149 



4yy HOME ; 'JTfl8 4 

410 INPUT '-NPlME OF M«P FILE^ 

420 OS = CHRS <4> 

430 PRINT OSj“OPEN MfiP” 

440 PRINT DSi “WRITE “;MS;".M«P“ 

450 PRINT L: PRINT W 
460 PRINT Uiz PRINT LV 
470 PRINT HL<0#0,0> 

460 FOR I = 1 TO L 
430 FOR J = I TO W 
500 PRINT HL< I^J,0> 

510 IF HLa,.J,0> = 0 THEN 550 
520 FOR K = 1 TO HL<Ii.J.0> 

530 PRINT HLa,J,k > 

540 NEXT K 
550 N0<T J 
t‘60 NEXT I 
570 PRINT OS; “CLOSE 

560 PRINT : PRINT : PRINT : PRINT DONE"; PRINT : PRINT 

590 INPUT “00 VOU W6NT TO CONTINUE INPUTTING OpTfl'?’ <V/N> '• iZS 

600 IF ZS = "V" THEN 200 

610 ijCiTO 30 

620 OS = CHRS (;4> 

630 PRINT OS; “OPEN ";MS;“.MPiP" 

640 PRINT OS;" READ ";MS;“.MPP" 

650 INPUT L: INPUT W 
660 INPUT LX: INPUT LV 
670 INPUT HL<0#0>0> 

630 FOR I = 1 TO L 
630 FOR J = 1 TO W 
700 INPUT HLa,J#0> 

710 IF HLa>J,0> = 0 THEN 750 
720 FOR K = 1 TO HL(I#J*0i'' 

730 INPUT HL< I^J,K> 

740 NE>iT K 

750 NEXT J 

760 NEKT I 

770 PRINT OS; “CLOSE “ ;MS; . NCiP" 

730 RETIJP.N 

790 PEN MHP CPTOL06 FILE READER 
300 HONE ; ‘■•JT H6 2 

610 PRINT TPB< li:;';“4*4- Mpp PEQDER 4*4" 

620 PRINT : PRINT 
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:.-:30. INPUT "NiaMe OP MfiP FILE? “ ;HS 
840 INPUT “PRINTER ON? <V/N> “;PS 
850 0« = CHRS <4> 

860 PRINT DSi“OPEN 

870 PRINT OSi“REflO “;H^;“,MAP“ 

880 INPUT L; INPUT W 
890 INPUT INPUT LV 

800 INPUT HL's0»0»0> 

810 FOR I = 1 TO L 

820 FOR J s 1 TO W 

830 INP*JT HL<I^J.0> 

840 IF HL<I#Ja 3) = 0 THEN 880 
850 FOR K s 1 TO Hta>J,0> 

860 INPUT HUI,.J,K> 

870 NEXT K 
880 NEXT J 
890 NEXT I 

1000 PRINT OS; “CLOSE “;HS;",HOP“ 

1010 HOME 

1020 IF as = "V“ THEN PR# I 
1030 PRINT “NflHE OF MPP: “;MS 

1040 PRU-<T “LONER LEFT ORID 30UPRE: “;LX;“ “ ;LV 
1050 PRINT “SIZE OF HPP; “:L;“ BV “;W 
1060 PRINT “8PSE ELEUPTION: ";HL<0.0^0> 

1070 PRINT 
1080 FOR I = I TO L 
1090 FOR J = 1 TO N 
1100 PRINT I>J,HL< I,J,0> 

1110 IF HLa,J,0> = 0 THEN 1160 
1120 FOR K s 1 TO HL<I^J,0::' 

1 130 PR I NT HL< I # J > i “ “ ; 

1140 ne*:t k 

1150 PRINT 

1 1 60 NEXT J 

1170 PRINT ; NEXT I 

1130 PR# 0 

1190 PRINT ; PRINT 

1200 INPUT “HIT RETURN TO CONTINUE” :XS 
I2l0 bOTO 30 
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E. HILL HAKER PROGRAM 



10 


REfi ■ 


HILL OPTP FILE HPKEP."'REPnER 


20 


DIH H<200#3),K';:3> 


30 K< 1 > s 5 


:K<2> = 3iK<3> » 14:K<4;;. = 13 


40 K<5> = 22;K(6> = 24;K<7> = 29:K<3> = 33 


50 


ftEH 


HILL OPTP FILE MPKER 


60 


HOHE : 


UTPB 2 


70 


PRINT " 


MENU: “ 


30 


PRINT " 


1 HILL MPKER" 


30 


PRINT "2 HILL REHOER" 


100 


PRINT 


"0 QUIT" 


110 


PRINT 


: INPUT "WHICH? ";X 


120 


IF X = 


» 0 THEN END 


130 


ON X I50TO 150,320 


140 


REH 


HILL OPTP FILE MPKER 


150 


HOHE ; 


: UTPB 2 


160 


PRINT 


TPBC hill MPKER 


170 


PRINT 


: PRINT 


130 


INPUT 


"NEW OR OLD FILE? 'rN.'O ’' " 


130 


IF PS 


= "N" THEN 230 



60SUB 560 

1‘10 INPUT "WHICH HILL 00 VOU WQNT tq ?T«RT WITH'? ";II 

220 GOTO 260 

230 WUME ; UTH8 4 

240 II a 1 



250 


INPUT 


"NUMBER OF HILL3? ";N 


260 


FOR I 


= II TO N 


270 


PRINT 


"ENTER THE OPTP FOR HILL ";I 


230 


INPUT 


"XCOQRO 13 ";H*;:i,l> 


290 


INPUT 


"VCOORD 13 ";H<I,2) 


300 


INPUT 


"PEPK ELEUPTION IS ";H<:i,3> 


310 


INPUT 


"ORIENTPTION PN6LE IS ":H''T.4-' 


320 


INPUT 


" ECCEHTR I C I TV 1 3 " ; H< 1 , 5 


7.30 


INPUT 


"SPREPD IS ";H<I,6;' 


340 


INPUT 


"MPX HEIGHT IS ";H<I,7:) 


350 


INPUT 


"CUT HEIGHT IS ";H<I>3> 


360 


PRINT 


; PRINT ; INPUT "SPUE'? CV.'N> ";2S 


370 


IF 2S 


= "V" THEN 390 


330 


NEXT i 


[ 


7-90 


HOHE ; 


: UTPB 4 


400 


OS = i: 


7HRS >;:4::' 


*10 


INPUT 


"NPME OF FILS'? ";NS 


420 


PRINT 


OS;"OPEH ";NS;“. HILLS, L34" 


430 


PRINT 


OS;" WRITE " ;NS;". HILLS, R0" 


440 


PRINT 


N 


+50 


FOR I 


= 1 TO N 


+60 


PRINT 


OS; "WRITE ", -NS; ".HILLS, P"il 
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470 FOR J s 1 TO 8 

480 PRINT 

490 NEXT J 

500 next I 

510 PRINT DSi“CLOSE “ ;N*;“ . HILLS" 

520 PRINT ; PRINT : PRINT DONE" 

525 PRINT : PRINT 

530 INP«JT "DO VOU MfiNT TO CONTTNOE INPUTTINe ORTR^? <V/N> ";As 
540 IF QS = "V" THEN 210 

550 GOTO 60 

560 HOME 

570 OS - CHRS <4> 



580 


INPUT "NRHE OF FILE^ ";N;S 




590 


PRINT OSj"OPEN " ;NS;“, HILLS. 


L34" 


600 


PRINT OS;"REfiD " ;N$; HILLS, 


R0" 


610 


INPUT N 




620 


FOR I = 1 TO N 




630 


PRINT OSi-REHD ";N«;“. HILLS, 


R" :T 


640 


FOR J = 1 TO 8 




650 


INPUT WI,J) 




660 


NEXT J 




670 


NEXT I 




530 


PRINT OS; "CLOSE ";NS;". HILLS 




690 


HOME 




P00 


PRINT "NUneER OF HILLS: " :N: 


PR IN 


710 


M = 1 




720 


FOR I = 1 TO N 




730 


IF M = 11 THEN PRINT :H = 1 




740 


PRINT I; 




750 


FOR J = 1 TO 8 




760 


PRINT TP6< K< J > >;H< I ,.J :■ ; 




770 


N0<T .J 




730 


PRINT 




790 


M = M -f 1 




800 


NEXT I 




810 


RETURN 
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:5ii0 ftEH HILL OftTP FILE REPDER 

330 RR = 0 

340 HOME ; UTPB 2 

350 PRINT TP8< HILL REPOEP 

360 PRINT ; PRINT 
370 O'S = CHRS •' 4 .:' 

360 INPUT “NPMe OF MPP'? “^-Nf 

360 INPUT "PRINTER ON*? <V/N> ";2'S 

300 IF Z£ - "V" THEN PR = I 

310 PRINT OS; “OPEN “;NS;". HILLS. L34" 

•320 PRINT OS;“REPD ";NS;“. HILLS. R0" 

330 INPUT N 

640 FOR I =* I TO N 

350 PRINT OS;“REPO " ;NS;“. HILLS. R";I 
:360 FOR J = 1 TO 3 
370 INPUT H<I..j;> 

660 NEXT J 
360 NEXT I 

1000 PRINT OS;"CLOSE " ;NS;“ . HILLS" 

1010 HOME 

1020 IF PR = 1 THEN PR# 1 

1030 PRINT ; PRINT NS;" HILL OPTP": PRINT 

1040 PRINT "NUMBER OF HILLS: ";N 

1050 PRINT 

1060 M = 1 

1070 FOR I = 1 TO N 

1080 IF M = 11 THEN PRINT :M = t 

1060 PRINT I; 

1100 FOR J s 1 TO 3 
1110 KK - K*r.J> + 3 * J 

1 120 IF PR = 1 THEN POKE 36. KK: PRINT hKT,,j-.;. 0OTO 1140 

1 1 30 PR I NT TPBC K V J > > ; HC I . -.1 ) i 

1 1 40 r4EXT J 

1150 PRINT 

1160 M = M * 1 

1170 NEXT I 

1 130 PR# 0 

1160 INPUT "HIT RETURN TO CONTINUE” ;XS 
1200 uOTO 6y 
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p. 300DS MAKER PROGRAM 



10 


REM 


WOODS CPTPLOG MPKER'- 


20 


DIM FL 


< 10, 10,9 ) 


30 


HOME : 


UTP8 3 


40 


PRINT 


"MENU: " 


50 


PRINT 


"1 WOODS MPKER" 


60 


PRINT 


"2 WOODS REPDER" 


?0 


PRINT 


"0 QUIT" 


30 


PRINT 


; INPUT "WHICH'? "iX 


:S0 


IF X = 


: 0 THEN END 


100 


ON X 


GOTO 120,770 


110 


REM 


WOODS CPTPLOG MPKER 


120 


HOME 


: UTPB 2 



130 PRINT Tfi8< 11 WOODS NftKER 

140 PRINT ; PRINT 

150 INPUT "NEW OR OLD MfiP*? <N.'-0 'J " :?-5 
160 IF ZS = "N*' THEN IT = 1:JJ = 1; SOTO 220 
170 INPtJT ’‘NRHE OF MOP? '';MS 
160 0$ = CHRS ':.4> 

1 90 SOSUB 600 

200 INPUT "ENTER STfiRTTNG GRID SOUfiRE a,J> ";II,JJ 
210 SOTO 240 

220 INPUT "SIZE OF MHP <1000 M ?QUPRES >'^ ":L..W 

230 INPUT "ENTER LOWER LEFT GRID SQUPRE CEG. 50,93) ";LK,LV 

240 FOR I = II TO L 

250 FOR J = JJ TO W 

260 PRINT ; PRINT 

270 PRINT "FOR GRID SQUPRE ":Ij",";J 

280 INPUT "NUMBER OF FORESTS'? " ;FL< I ,J,0 ):KK = FL<I,J,0) 

290 IF KK = 0 THEN 370 

300 PRINT "ENTER THE ";KK;" FOPFSTS FOR GRID SOUQRE 
310 FOR K = 1 TO KK 

320 PRINT "FOREST ";K;" IS INPUT FL< I ,.J,K :• 

330 NEXT !<■ 

340 JJ = I 

350 PRINT : INPUT "SPUE?’ < V.-'N)" ;p$ 

360 IF PS = "V" THEN 390 
370 NEXT J 
330 t€XT I 
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1.30 HOME : UTP8 4 

400 INPUT “NPME OF MfiP FILE? “ ;M« 

410 OS s CHRS <4> 

420 PRINT OSi“OPEN " ;HSi “.WOODS'* 

430 PRINT OSj“WRITE “;MS;“. WOODS- 
440 PRINT L; PRINT W 
450 PRINT LX: PRINT L’l* 

460 FOR I = 1 TO L 

470 FOR J s 1 TO W 

460 PRINT FL< I,J,0> 

490 IF FL<I>J,0> = 0 THEN 530 
500 FOR K = 1 TO FL<I^J,0> 

510 PRINT FL<I,J,K> 

520 NEXT k; 

530 NEXT J 
540 NEXT I 

550 PRINT os; “CLOSE “ ;MS;“. WOODS" 

560 PRINT : PRINT : PRINT : PRINT "... DONE'*: PRINT : PRINT 

570 INPUT “DO 't'OLI WPNT TO CONTINUE TNPUTTIMS OPTO-^ < 'I'.'N :■ “ ;2S 

580 IF ZS = “V“ THEN 20O 

590 GOTO 30 

600 OS = CHRS <4> 

610 PRINT OSi“OPEN “;MS;“. WOODS" 

620 PRINT OS;“REPD “;MS;“. WOODS'* 

630 INPUT L; INPUT W 
640 INPUT LX: INPUT LV 
650 FOR I = 1 TO L 
660 FOR J = 1 TO W 
670 INPUT FL(I»J^0> 

660 IF FLa,J,0::> = 0 THEN 720 
690 FOR K = 1 TO FL’kI»J^0.> 

700 INPUT FLa^J^K) 

’10 next K 

720 NEXT J 
T30 next I 

’40 print DSj “CLOSE “;MS;“, WOODS" 

750 RETURN 
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'sei «EH WOODS CWTfiLOG FILS ftSPOER 
770 HOME ; UTfie 2 

r30 PRINT Tfiec WOODS RECiDEP 

790 PRINT ; PRINT 
900 INPUT "NRME OF MfiP FILE? " :HS 
910 INPUT "PRINTER ON'?' CV/n;:- ";R:r 
920 OS - CHRS <4> 

930 PRINT OS; "OPEN ";HS;‘‘. WOODS" 

940 PRINT OS;"REfiO " ;MS: " .WOODS" 

950 INPUT L; INPUT W 

9S0 INPUT LK: INPUT LV 
370 FOR I = I TO L 
990 FOR J = 1 TO W 
990 INPUT FL<I,J-0> 

900 IF FL<I^J,0> = 0 THEN :940 
910 FOR K = 1 TO FLU,J.0> 

920 INPUT FL<I,J,K> 

930 ne:<t k 

940 NEXT J 

950 NEXT I 

:J60 PRINT OS; "CLOSE " :MS;". WOODS" 

970 HOME 

980 IF AS = ”V“ THEN PR# I 
990 PRINT “NAME OF MAP: ";MS 

1000 PRINT "LOWER LEFT GRID SOUAREr “;LX.;" " ;l.v 

1010 PRINT "SIZE OF MAP: ";L;" BV " ;W 

1020 PRINT 

1030 FOR I = 1 TO L 

1040 FOR J = 1 TO W 

1 050 PRINT I , J , FLC I , J ,0 > 

1090 IF FL<I^J*0> = 0 THEN 1110 
1070 FOR K = I TO FL':;I*J#0> 

1 080 PR I NT FL( I , J ,K > ; " " ; 

1090 NEXT K 

H00 PRINT 

1110 NEXT J 

1120 PRINT : NEXT I 

1130 PR# 0 

1140 PRINT ; PRINT 

1150 INPUT "HIT RETURN TO CONTINUE" :X? 

1190 GOTO 30 
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G 



FOREST MAKER BROGRAM 



lti REH FOREST MRKER/RERDER 
20 DIM F< 120*8 )»K*' S ) 

30 K<1> = 5:K<2> = 12:K<3> = 20;KX4 > = 2R 
40 Kv5> = 30:K<6> = 38 
50 HOME : UTR8 2 
80 PRINT "MENU: " 

T0 PRINT "1 FOREST M(:iKER“ 

80 PRINT "2 FOREST READER" 

S0 PRINT “0 QUIT" 

100 PRINT ; INPUT "WHICH? ";X 

110 IF X = 0 THEN END 

120 ON X 80T0 140,300 

130 REM FOREST DATA FILE MAKER 

140 HOME : UTA8 2 

1 50 PRINT TA8< 10>;"4** FOREST MAKER 
180 PRINT : PRINT 

1 70 INPUT "NEW OR OLD FILE? ''N/0> ";AS 
130 IF As = "N" THEN 220 
190 80SU8 530 

£00 INPUT "WHICH FOREST DO VOU WANT TO 3TART WITH 

210 80T0 250 

220 HCiHE ; UTA8 4 

230 11=1 



240 


INPUT 


"NUMBER OF FORESTS'? " ;N 


250 


FOR I 


= II TO N 


280 


PRINT 


"ENTER THE DATA FOR FOREST ";T 


270 


INPUT 


"XCOORD IS "*F<I*n 


280 


INPUT 


"VCOORD IS ";F-ri,2> 


290 


INPUT 


"HEI8HT OF TREES IS ";FvI*3::' 


300 


INPUT 


"ORIENTATION AN6J.E IS ";F<'I.4> 


310 


INPUT 


"SEMI-MAJOR AXIS LENGTH IS " ;F*:; 1 ,5 


•320 


INPUT 


"SEN I -MI NOR AXIS LENGTH IS " :F< I .8 


330 


PRINT 


: PRINT : INPUT "SAUE'? -: V--N> ";2S 


340 


IF 25 


= "V" THEN 360 


350 


NEXT 


I 


360 


HOME 


; UTA8 4 


370 


OS = 


CHRS ■:'4 ;i 
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38€i INPUT “NfiMft OF FILE? “ ;NS 

390 PRINT DSi"OPEN NS; “.FORESTS. L3S“ 

4.00 PRINT 0S;“ WRITE NS; ".FORESTS, R0“ 

410 PRINT H 

4i0 FOR I = I TO N 

430 PRINT OS;“WRTTE NS; “.FORESTS, R“;T 
440 FOR J = I TO S 
450 PRINT F<I,J> 

400 next J 
470 N£:<T I 

490 PRINT OS .-“CLOSE NS; “ .FORESTS" 

490 PRINT ; PRINT : PRINT "... DONE" 

495 PRINT : PRINT 

500 INPUT "00 VOU WftNT TO CONTINUE INPUTTING nWTP'? < V/N > 
510 IF QS s "V" THEN 2O0 
520 GOTO 50 
530 HONE 

540 OS = CHRS <4 ;. 

550 INPUT “NfiHE OF FILS'? ";NS 

560 PRINT OS; “OPEN NS; “.FORESTS, 139" 

570 PRINT OS;"RERD "; NS; ".FORESTS, R0“ 

530 INPUT N 



590 FOR I = 1 TO N 

600 PRINT DS;“RE0O “; NS; ".FORESTS, R":I 
610 FOR J = 1 TO 6 
620 INPUT F<I,J> 

630 NEXT J 
640 N0<T I 



650 PRINT OS; “CLOSE " ;NS;“ . FORESTS" 



660 HOME 

670 PRINT NS;“ FORESTS LIST": PRINT 



680 PRINT “NUMeER OF FORESTS: ";N: PRINT 
690 N = I 



700 FOR I = 1 TO N 

710 IF M = 11 THEN PRINT :M = I. 

720 PRINT I; 

730 FOR J = 1 TO 6 

740 PRINT T0B< K< J>>;F< I,J >; 

750 NEXT J 

760 PRINT 

770 M = M + 1 

730 NEXT I 

790 PET URN 



;0S 
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80V3 REM FOREST DOTP FILE READER 
310 HOME ; UTAB 2 
320 PR s 0 

330 PRINT TA8< 10 FOREST READER 

340 PRINT ; PRINT 
350 O'-f * CHR* <4> 

360 INPUT “NAME OF MAP? “;N$ 

370 INPUT “PRINTER ON? <V/N> “;Z$ 

3:30 IF ZS = “V THEN PR = 1 

:390 PRINT DS;“OPEN ";N$i“. FORESTS. L3S" 

300 PRINT OS; “READ ";NS;'‘. FORESTS, R0“ 

310 INPUT N 

320 FOR I s 1 TO N 

330 PRINT OS;“REAO NS; “.FORESTS. R“;I 

;340 FOR J s 1 TO 6 
350 INPUT F<I,J> 

360 NEXT J 
370 NEXT I 

380 PRINT DS; “CLOSE “; NS FORESTS “ 

390 HOME 

1000 IF PR s 1 THEN PR# I 

1010 PRINT : PRINT NS; “ FOREST DATA"; PRINT 

1020 PRINT “NUMBER OF FORESTS ";N 

1030 PRINT 

1040 M s 1 

1050 FOR I = I TO N 

1060 IF M * 11 THEN PRINT tM = t 

1070 PRINT I; 

1080 FOR J s 1 TO S 
1030 KK s KXJ> + 3 * J 

1100 IF PR a I THEN POKE 36, KK: PRINT F-' I.. SOTO 1120 
1 1 10 PRINT TAB< K< J > > ;F< I , J > ; 

1120 NEXT J 
1130 PRINT 
1140 M a M + 1 
1 150 ^€XT I 
1160 PR# 0 

1170 INPUT “HIT RETURN TO CONTINUE” ;XF 
1 130 30T0 50 
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H. TSRHAIH MODEL PROGRAM 



10 


REM 


UPDATED 02 


PE8 82 




20 


REH 


D0MSTfiC TERRRIH MODEL 




30 


DIM 


XB< 5 >/t'8< 5 > 


,XR< 5 > 5 5 


, ?. 0 •' . 8 V< 5 * 1 0 > . RX< 5*10 ;• 


|35 


DIM 


ZB<5>.2R<5> 


,HL< 10* 10, 12 ■> 




40 


DIM 


8P<5>,RP<5-' 


,S8<5>*3R(5>,ftV<'5 


, 1 0 


50 


DIM 


eF<5>,RF<5> 


*eT<5*5>*RT<5,5)*^ 


DB<5')*0R<5> 


60 


DIM 


8R< 5 >,RR< 5 > 






70 


DIM 


BK5:>,R1<5;:' 







T = 0!OT = 30! HT » 1500 
90 OS = CHRS <4> 

100 P = 0.99:Q = 0.93 
110 PI = 3. 141592654; PR = 0 
120 0 » 2. 4; 8 = 0.03 
130 BE = 0.3: RE = 0.3 
140 8M = .5:RH = .5 
150 K.K = 0 



160 


NK = MT 


/ OT 


165 


22 = 1 




170 


DIM 80< NK > *R0< NK ■> 


171 


HOME 




172 


INPUT 


"NHME OF MOP FILE'^ “ ;MS 


173 


PRINT 


OS; “OPEN “;MS;“.MOP“ 


174 


PRINT 


OS * “REPO “ ; MS ; . MPP “ 


175 


INPUT 


L: INPUT W 


176 


INPUT 


LX: I NPt JT Ut> 


177 


INPUT 


HL*: 0*0*0 > 


178 


FOR I 


= 1 TO L 


179 


FOR J 


= 1 TO W 


180 


INPUT 


HL< I*J*0> 


131 


IF HL< 


I*J.0;' = 0 THEN 184 


182 


FOR K 


= 1 TO HL< T*J*0:;' 


183 


INPUT 


HLt; i*.j*K.:;' 


184 


NEXT ; 


i NEXT : NEXT 


185 


PRINT 


OS; “CLOSE “;MS;“.HPP“ 


136 


LX = l:= 


< - 1:LV = LV - 1 


187 


PRINT 


OS; “OPEN “; MS; ".HILLS* L34 


138 


PRINT 


OS:“REPD “ ;MS;“. HILLS. RP" 


l39 


INPUT 


NH: OIM H<NH*8:) 


190 


FOR I 


= 1 TO NH 


191 


PRINT 


OS;“REPO “;MS;“. HILLS* R“ ; 


192 


FOR J 


= 1 TO 8 


193 


INPUT 


H( I *.J > 


194 


NEXT 


: NEXT 


195 


PRINT 


OS CLOSE “ ; MS ; “ . H T L L S " 
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I 8a HOME 

i89 Input “NfiME OF FILE FOR SLUE OPTQ*- “ ;HS 
200 PRINT OS; "OPEN NS; “ .BLUE” 

210 PRINT OS;“RERO ";NS;“, BLUE" 

220 INPUT M 

230 FOR I = 1 TO M 

240 INPUT B1 

250 FOR J s I TO B1 

260 INPUT BX<I,J>: INHJT BV< I > 

270 NEXT ; NEXT 

280 PRINT OS; "CLOSE “;NS;".BLUE’‘ 

280 FOR I s 1 TO M:XB<t> = BX< 1 ,1 >: VB< D = BVa.DiBPM) = t. 

300 PRINT 81: FOR JJ - I TO 81: PRINT BXCI ,J..i >-8V< I ,JJ >: NEXT JJ 

310 PRINT "SIZE OF BLUE UNIT ”;I:r INPUT SBC T 
320 eu I > = SBCI > 

330 B0< 0 > = B0< 0 > + SBC I > 

340 BR< I > = 0 

341 X s XBC I>:V = vec I > 

.7.42 X4 s I NT <XB<I> 1000) - i.X:V4 = TNT 'VBCr> 1000) - LV 

343 603U8 I960 

344 ZSC I ;:■ = 2 

345 PRINT "ELEWION OF SLUE IS " ;Z 

346 FOR II = 1 TO HLCX4,V4,0>:H2 = HL<X4pV4,U );HCH2,0) = 0: NEXT 
350 NEXT 

360 PRINT 

370 INPUT "NAME OF FILE FOR RED OpTQ? •• ;NS 
380 PRINT OS; "OPEN ";NS;“. RED" 

380 PRINT OS;"RERO "; NS; ".RED" 

400 INPUT N 

*10 FOR I = 1 TO N 

420 INPUT Rl 

*30 FOR J = 1 TO Rl 

440 INPUT RXCI.J): INPUT P.VC I ,.J 

*50 NEXT : NEXT 

*60 PRINT OS; "CLOSE ";NS;". RED" 

*70 FOR I s 1 TO N:XR<n = RXC I , I "'; VRC I = RV< I . f- : RPC I > = 1 
480 PRINT "SIZE OF RED UNIT ";I;: INPUT SRC I .• 

*80 Ric I > = 3R< I 

500 R0C. 0 > = R0C 0 > + SRC I ) 

510 RRC D = 5 



511 : 


X = XRC 


I ;:';V 


= VRC I 


’’t 










512 


X4 = I 


NT c: 


XRC I > 


1000 :i - 


LX:V4 = 


I NT > VRCr> 


.•■■■ 1. 000 ) 


- LV 


513 


803UB 


i;360 














514 


ZRCI> = 
















51d 


PRINT 


"ELE' 


URTION 


OF RED " 


;I?" IS " 








4 516 


FOR 11 


= 1 


TO HLC 


X4,V4,0 >1 


tH2 = HLC 


X*,V4.I1 ):HC 


H2 » 0 •• ~ 


0: NEXT 


520 


NEXT I 
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■530 PRINT : INPUT “00 VOU NpNT TP*= PPTMTPR <- v./n > “ ;Ps 

535 IF PS = “V* THEN PR = 1 

540 HCiME : UTP6 3 

550 IF PR * 1 THEN PR# I 

560 PRINT 

570 PRINT ”TIME“i: POKE -36.7? PRTNT “5UJE'‘^; °0KE -36,20: PRINT “RED“ 
560 PR# 0 

590 IF Z2 * 1 THEN i7»HSUB t760: GOTO 610 

600 G0SU6 1490 

610 REH MPIN PR06RPH 

620 T = T + OT 

630 IF r < s NT THEN 640 

635 PRINT “ENO OF 8PTTLE DUE TO TINF"* G06UB 1540? SNO 

640 G0SU6 790 

650 G*3SU8 1170 

660 KK = KK + 1 

670 80<KK> = 0:R0<KK') = 0 

660 FOR I — 1 TO M 

690 60< KK > = 80< KK > SB< I > 

700 NEXT I 

"10 FOR J s I TO N 

720 R0(Kk> = R0<:KK> + SR<J> 

730 NEXT J 

740 IF 22 = 1 THEN GQSUB 1760: GOTO 
750 I50SU6 1490 

760 IF 80<KK> > BE ^ 80<0 ) THEN 770 

765 PRINT : PRINT “ENO OF BPTTLE DUE TO PTTRITION OF BLUE" 

766 iSOSUB 1540: ENO 

770 IF R0<KK> > RE 4 R0<0> THEN 790 

775 PRINT : PRINT “END OF RPTTLE DUE I’D PTTRITION OP RED" 

778 GOSUe 1540: END 
730 GOTO 620 

790 REN NOUENENT SUBROUTINE 
300 FOR I = 1 TO N 
310 J = RPC I > 

320 OV = RVC I , J + r> - VRC I ■> 

330 OX = RXt I , J 1 '■.< — XRC I .> 

:340 IF OX = 0 PNO OV > 0 THEN QN = O’! pt Gt''TO 996i 

350 IF OX = 0 PNO OV < 0 THEN PN = 3 * PI 2? GOTO 390 

360 PN = PTN CDV / riX> 

370 IF OX < 0 THEN PN = PI + PN 
380 01 = OT 4 RR< I :> * COS CPN > 

390 02 = OT 4 RRt, I> ^ SIN CPN;:' 

900 OH = 30R CDl 2 4 02 •••' 2"' 

910 01 = 9QF; ( OX 2 4 OV 2 

920 IF 01 > = OH THEN 950 

930 01 = OX: 02 = DV 
340 RPC I > s J 4 1 
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10£i»5 ) - LV 



850 VR<I> = IMT CVP:<;i> + 02 

860 XR< I :• = TNT ( KR< I > 01 ) 

861 X = XR< I >:V = VR< I > 

862 X4 = I NT <XR(I> •••' 1000> - U<:V4 = I NT •: VR< T > 

;963 60SUB I860 

:864 2R< I = 2 

865 PRINT "ELEUPTION OF RED "il;*' TS " ;2 

:366 FOR II = 1 TO HL< X4,V4,0 ):H2 = l-H.'::X4,V4,Il >:WH2,0) = 0: NEXT 
870 NEXT 

880 FOR I » 1 TO M 
890 J s ep< i> 

1000 Dv = 8 v< i,j + 1 > - ve< T > 

1010 OX = BX< I,J + 1 > - XB<’ I 

1020 IF OX = 0 PND DV > 0 THEN PN = PT / 2: 80T0 1069 
1030 IF DX s 0 OND DV < 0 THEN PN s 3 4 PI / 2: 60T0 I860 
1040 PN = PTN <DV / DX> 

1050 IF OX < 0 THEN PN = PI + PN 
1060 01 = OT * BR<I> 4 COS <PN> 

1070 02 = OT 4 BR<I > 4 SIN <PN;;' 



02 

DV 



1080 OH = SQR <D1 
1080 01 = 30R COX 
1100 IF 01 > = OH THEN 1130 

1110 01 = DX:02 = OV 



^24 

-s 2 + 



2 ■' 



1120 6P< I 



= J 4 I 



1130 


VBC I) = 


INT CVBCn 4 02 ) 




1.140 


XB< I > = 


TNT CXBC n 4 01 ) 




1150 


NEXT 






1160 


RETURN 






1170 


PEH PTTRITTON SUBROUTINE 




1180 


FOR I = 


1 TO M:8FCI> = 0;OBCI> = O; 


NEXT 


1 180 


FOR J = 


1 TO N:RFCJ> = 0:ORCJ> - 0: 


NEXT 


1200 


FOR I = 


1 TO M 




1210 


FOR J = 


1 TO N 





1220 BTCI^J) = 0:RTCJ,n = 8 
1 230 RG = SOR C C XSC I ) “ XR< J > > 



2 4 .r VBC I > - VRC J > :■ 



1240 


IF RG 


> 2700 THEN 


1310 








1250 RR( I 


w 










1251 


REH 


CHECK LINE OF 


SIGHT 








1252 


XI = 


TNT <X8<I> 


1000 ;■ - LX;Vt = 


TNT 


■' VB< T > ••• 


1000O - LV 


1253 X2 = 


INT CXRCJ> 


1000> - LX:V2 = 


INT 


CVRC.J-;- 


1000 > - LV 



1254 

1255 

1256 K4 = 



3 = X2 - X1:V3 = V2 - vt 
FOR II =0 TO 



PBS CX3> 

II 4 II SON C X3 !•' 

1257 FOR J1 =0 TO PBS 'rV3> 

1253 V4 = VI 4 Ji * SON CV3’:- 

1261 PRINT “CHECKING GRID SOUPRE " ;X4;“ , “ :V*i ; PHR BLUE RED 
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12S2 


6osue 


I960 












1263 


N0<T ; 


NEXT 












1264 


FOR 11 


s 1 TO NH:H< 11. 


>0J 


' = 0; 1 


NEXT 


1265 


REM 1 


LOS IS 


OK 










1270 


BF< I > = 


8F< I > 


1 










1280 


RF<J> = 


RF<J> 


+ 1 










1290 


8T( I,J> 


- M ^ 


<i - 


< P -• 




3B< I > > 


3R^ -.1 ' 


1300 


RTc-J,I> 


= B # 


<1 - 


0 > 




SRiC J > ) 


# S6< I > 


1310 


NEI^r : 


NEXT 












1320 


FOR I 


= 1 TO 


M 










13-30 


FOR J 


- I TO 


N 










1340 


IF RF< 


•J> = 0 


OR BF( I > 


S 


0 THEN 


1370 


1350 


DBCI) = 


DB( I > 


+ BT-:; 


I,J 


y / 


' rf<j:^ 




1 360 


OR<J> = 


i OR<.J> 


■K 9J< 


JO 


> y 


•' BFvi ;- 




1370 


NEXT : 


NEXT 












1380 


FOR I 


s 1 TO 


M 










1390 


SBf I > = 


^ 3Ba> 


- ORC. 


n 








1400 


IF S8< 


I > < 0 


THEN 


se< 


I'.' 


= 0 




1410 


IF S6< 


i:> <• = 


= BM 


• 81 


< I 


> THEN 


BR( I > = 8 


1420 


NEXT 














14-30 


FOR J 


- 1 TO 


N 










1440 


SR<J> = 


* 3R<J> 


- 0R( 


J> 








1450 


IF 3R< 


J) < 0 


THEN 


3R< 


J> 


= 0 ' 




1460 


IF 3R< 


I> < ^ 


= RM 4 R1 


< I 


:> THEN 


RR< I ■) = 


1470 


NEXT 















14Si£| return 

1490 REM PRINTOUT/SUMHfiRV 
1500 IF PR - 1 THEN PR# I 

1510 PRINT T;: POKE PRINT B0<KK>;: POKE 36^^20: PRINT R0<KK > 

1520 PR# 0 
1530 RETURN 

1340 REM 3HUE RESULTS IN H FILE 

1550 PRINT : INPUT “DO VOU WONT R FILE OF THIS RUN^S ORTR? ^V/N> “jRS 
1500 IF RS = "H" THEN 1750 
1570 OS = CHRS <4> 

15S0 INPUT “NRHE OF FILE FOR FORCE LEUP. ORTR-? '• ;FS 
1590 PRINT OS; “OPEN " ;FS;".COEFS" 

1600 PRINT DS;“WRITE '* ;FS;*m:OEFR" 

1610 PRINT R: PRINT 8: PRINT P*. PRINT 0 
1620 PRINT OS; “CLOSE ” ;FS;“ . COEFS“ 

16-30 PRINT OS; “OPEN “ ;FS; RESULTS" 

1840 PRINT OS; "WRITE “ ;FS; RESULTS" 

16-50 PRINT NK + 1; PRINT HK + 1 

i860 FOR I - 0 TO NK 

1670 T = I » 3 

1690 PRINT T: PRINT 800 

1690 NEI'iT I 
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1700 FOR I s 0 TO Nk' 

1710 T = I * 3 

1720 PRINT T: PRINT R0< I > 

1730 NEKT I 

1740 PRINT OS,-“CLOSE “ RESULTS’* 

1750 RETURN 

1760 REM PRINTOUT/SUMMPRV 

1770 HOME : UTfiB 3 

1730 IF PR = 1 THEN PR# 1 

1790 PRINT “TIME IS “;T 

1300 PRINT : PRINT 

1310 PRINT “BLUE FORi:ES“ 

1320 PRINT “UNIT”; TPB< 7>;“XCOORD“; TCiB< i5>;“'yC00R0“ ; TP8< 23 );“Si: 
1330 FOR I = 1 TO M 

1340 PRINT “ “;I; Tfte< 7:);X8< I >; TR8< 15 >;‘r*8< ! >; TP8( 24>;SB< I -* 

1350 NEXT 

I860 PRINT : PRINT “RED FORCES” 

1370 PRINT “UNIT"; TP8< 7::-;“XC00R0"; TP8< 15:j;“‘t>C 00R0"; TP8< 23>;“Si: 
1380 FOR I = 1. TO N 

1330 PRINT “ “;I; TPB< 7>:XR<I>; TP8< 15:>;*yR<I>; TPB< 24>;SR<I> 

1800 NEXT 

1910 PRINT : PRINT 

1920 PR# 0 

1930 INPUT "HIT RETURN TO CONTINUE" ;Z? 

1940 RETURN 

1950 REM ELEUPTION COMPUTPTION 

I960 CR ~ - 99999999 

1970 2 = HL<0ii0>0:) 

I960 FOR HI s 1 TO HL<X4,V4,0 :- 

1990 H2 = HL(X4,V4,Hn 

2000 IF H<H2,0) = 1 THEN 2150 

2010 X3 = X - H<H2,1> ^ 100:*T»S * *J - H<H2.2> 100 

2020 PI s LOO <H<H2^7‘;' < H<' H2,7 - 30 >> 

2030 SI = PI 4 <H<H2^5>> 2 

2040 Cl = H<H2,4> 4 PI 180 
2050 C2 = H<H2^8> ^ 2 

-060 Pi “ “CPl 4C cOS Cl ' 2 ■** 8 1 4 ^ SIN Cl ' ■*“ 2.-' 

2070 P2 — — <P1 4 i' SIN >' Cl > ) ''2 + 81 4 COS- '"Cl '■ > 2) •••' C2 

2030 P3 ~ <2 4 COS ';C1> 4 SIN *'01. •' 4 f'8l - pi. > C2 
2090 IF H< H2..7') > H<H2,3> THEN CR = L03 H<' H2 ^ 7 - H< H2,3V;. H-:' H2 

2100 01 = Pi 4 X3 .-s 2 4 P2 4 VS 2 + P3 4 VS + VS 
2110 H<H2,0 :' = 1 
2120 IF 01 <• CR THEN 2150 

2130 FI = H<H2,3> + H< H2,7> 4 £XP '::or:' - 1 > 

2140 IF FI > 2 THEN 2 = FT 
2150 NEXT 
2160 RETURN 
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